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XLVIT. The Effect of Small Amounts of Cold-Work on Y oung’s Modulus 
of Copper 


By A. D. N. Smrra 
Royal Aircraft Establishment, Farnborough* 


[Received December 5, 1952; Revised January 5, 1953] 


SUMMARY 


The effect on Young’s Modulus of copper of small amounts of plastic 
extension (1-8%) and subsequent low temperature annealing (50—-100°c) 
has been studied. The results confirm previous work that there is a 
reduction of about 11% in the modulus, the greater part of which is 
recoverable after prolonged annealing at 100°c. The reduction is attri- 
buted to the elastic oscillation of the relatively mobile screw components 
of dislocation loops, recovery being caused by the annihilation of such 
components in pairs by cross-slip, as suggested by Mott (1952). 

In an analysis of the recovery process it is assumed that the activation 
energy required to move the two screw components together is not unique 
but is distributed in a Gaussian manner about a mean. This mean value 
is found to be 25 000 cals/mole for O.F.H.C. copper extended by 1% and 
21 500 cals/mole for the same material extended by 8-4%. 


§ 1. INTRODUCTION 


As long ago as 1907 Beilby noticed that the elastic moduli of metals are 
reduced by small amounts of cold-work, although partial recovery may be 
effected by low temperature annealing, and this phenomenon has since 
been confirmed by many investigators (e.g. Kuntze 1928, Kawai 1930, 
Koster 1940). It has also been noted that while a small amount of 
deformation reduces the modulus, continuing cold-work may reverse this 
trend in some cases, and may actually cause the modulus to increase above 
its original value. While it seems clear that the latter effect is due to the 
development of preferred crystallographic orientation, coupled with the 
elastic anisotropy inherent in most metals, the reason for the former has 
not yet been established although three theories have been put forward, 
which will be examined in detail in § 4 below. Some experiments have 
been made to assess which, if any, of these three theories is correct. In 
these experiments a number of polycrystalline copper bars were stretched 
plastically by small amounts, Young’s Modulus, #, being measured 
dynamically before and after extension. The bars were then annealed for 
various times in the temperature range 50-100°c, the modulus being 
measured at sufficiently frequent intervals during the annealing period to 


establish a series of recovery curves. 
ee ES ee 
* Communicated by the Author. 
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§ 2. EXPERIMENTAL DETAILS 


Material : The copper used in these experiments was principally some 
oxygen-free, high-conductivity copper supplied by Messrs. Thomas 
Bolton Ltd., in the form of rolled sheet of cross section 6 in. x +in., from 
which specimens of dimensions approximately 6 in. x } in. x $ in. were cut. 
A typical maker’s analysis of material supplied to the same specification 
as that used here.is as follows: Cu 99-984%, Ag 0-002%, Fe 0-0028%, 
Ni 0-:003%, Pb 0-001%, Bi 0-0004%, O» 0:0045%. No trace was found of 
Se, Te, As or Sb. 

Three experiments were also made using spectrographically pure copper 
in the form of 6 mm diameter bars, supplied by Messrs. Johnson Matthey 
and Co. Ltd. 

Apparatus: Young’s modulus, #, was measured dynamically by 
determining the fundamental frequency of vibration in the ‘ free-free 
mode, using the relation 


0:0787 f2 
ii as 


where f is the natural frequency in cycles per second, / the length of the 
bar, p the density and k the radius of gyration of the section. 

The apparatus used to determine the natural frequency followed the 
design of Férster (1937). Briefly, the specimen was suspended horizontally 
by two fine cotton threads attached to the bar close to the nodes of the 
‘ free-free ’ mode, the other ends being attached to the coil of a moving-coil 
loudspeaker and to a gramophone pick-up, respectively. The loudspeaker 
which was driven from an accurately calibrated decade oscillator, into! 
vibrations in the bar which were detected by the pick-up, amplified and fed 
to a cathede ray oscilloscope. The oscillator was tuned so that the 
vibrations increased to a maximum, and the resonant frequency was then 
read directly from the decade dials. : 

Experimental Procedure : Kach specimen was first annealed in vacuo for 
3 hours at 500°c, after which treatment x-ray back reflexion photographs 
showed that recrystallization was complete, while the grain size eo one 
sample bar was about 0:06mm. The moduli of the annealed bars were 
measured and then each bar was given a small amount of plastic strain in a 
tensile testing machine. Of the O.F.H.C. bars, seven were niente b 
1° and seven by 8-4°%, while all the three spec. pure bars vere ee 4 
1%. The ends of the test piece were then sawn off to remove the ty. 
made by the grips of the machine. The modulus was measured a i a d 
then each bar was annealed in the range 47—100°c, the me a ee 
measured at suitable intervals during heat treatment. To ensure a cl “ 
control on the annealing temperature for the longer anneals, the bars it 
first heated up quickly by immersion in hot water at the . ro viet 
temperature, then transferred to a vacuum flask full of oil ee ihe i 
temperature, which was itself placed in a thermostatically-controlled cae 
The oil served the dual purpose of smoothing out the tamporature ae 
ations of the oven and reducing the oxidation of the Ae For pei 
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anneals (< 3 minutes), the specimen was kept in hot water for the whole of 
the annealing period. 

Finally when sufficient points on the recovery curve had been obtained, 
all the bars were given a prolonged anneal at about 100°c, after which 
treatment it was assumed that all had recovered to approximately the 
same degree. 

§ 3. RESULTS 


O.F.H.C. Copper, 1%, Extension: The recovery curves are shown in 
fig. 1, in which the dimensionless parameter (H,—E)/E, is plotted against 
log, ot, where H, was the final value of the modulus after a prolonged anneal 
of 70 hours at 105°c, and H was the modulus after ¢ hours at the given 
temperature. The error in this parameter due to uncertainty of measure- 
ment is probably about -+-0-003. 


Fig. 1 


O-0! ol | 10 100 HOURS 


ANNEALING TEMP. 


-2:0 -1-0 ° ie) 2:0 


LOG,, t 
(t MEASURED IN HOURS) 
Recovery curves for bars extended 1%: O.F.H.C. copper. 


It will be seen from fig. 1 that the recovery curve for the 52°c anneal 
is incomplete, owing to the accidental dropping of the specimen after it had 
been annealed for a few hours. It has been assumed however that the 
value of #, was not materially altered by this accident. 

O.F.H.C. Copper, 8:-4% Extension: The corresponding curves for the 
seven. bars in this group are shown in fig. 2, the final common anneal being 
for 230 hours at 100°c. 

Spectrographically Pure Copper, 1° Haxtension: Figure 3 shows the 
recovery curves for the three bars in this group, the final common anneal 
being for 70 hours at 105°c. 

Effect of Stress: One experiment was made to discover whether the 
recovery process was speeded up by the application of stress. Two bars 
were extended by 1%, and after measuring the modulus of each, one was 
driven at resonance by the loudspeaker, while the other was not, both bars 


2H2 
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being kept at the same temperature (18°c). The amplitude of vibration 
used in this experiment was about 4 times greater than that normally 
used to measure the modulus. After about five hours of this treatment 
the moduli of both bars were measured in the normal manner, and it was 
found that both had increased by about 0-8%. It was also found that 
increasing the vibration amplitude by a factor of four did not alter the 
modulus detectably. Hence the reduction in modulus due to cold work is 
not influenced by the magnitude of the stress nor by the length of time for 
which it is applied, under the conditions of this experiment. 


Fig. 2 
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Recovery curves for bars extended 8-4% : O.F.H.C. copper. 


Effect of Removal of Ends of Test-Piece : Since it was possible that the 
working introduced by sawing off the ends of the bars could have a marked 
effect on the modulus, the ends were sawn off an es, bar; it was 
found that # was reduced by 1° after one end was removed and by a 
further 0-5% after the other end was taken off. Thus the éffect was bee h 
smaller than the reduction of about 1 15% due to plastic axteneres “ 

Effect of Reduction in Temperature: One bar was cooled down to sub 
normal temperature after it had been worked, to see whether there ae 
substantial increase of the modulus associated with reduction in snipe 


ature. The amount of plastic extension was 1°, and the values of the 
modulus were : j | 


Fully annealed 12 700 kg/mm? 
Freshly worked at 19°C 11 440 ke/mm? 


—15°o 11 650 kg/mm? 

—41°o 11 840 ke/mm2 

After test was finished 19°c 11 530 mane 
g) 


poked SEE a Sa ! 12a 100 HOURS 


(ml Ope ee ee 


Ep 
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Thus, taking the mean of the values before and after cooling as the norm, a 
reduction of 60°C increased the modulus by about 3:1%. This are 
should be compared with the increase to be expected in the modulus of 
annealed copper over the same temperature range, given. by Koster (1948) 
as about 1:9%. Hence it appears that there is some increase of the 
modulus on cooling by 60°C below room temperature, but that it is only 
about 10°% of the original loss. 


ANNEALING TEMP. 
oor 
79°5 

52 


KEYS X 
1c} 
+ 


cate) 


LOG, t 
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Recovery curves for spec. pure copper extended 1%. 


§ 4. DiscussION 


Possible Causes of the Reduction in Modulus: Possible causes of the 
reduction and the subsequent recovery after low temperature annealing 
have been put forward by Zener (1948), Eshelby (1947) and Mott (1952). 

Zener attributes the reduction, together with the Bauschinger effect, to 
the viscous nature of freshly formed slip-bands. If this is the case the 
effect should be very sensitive to temperature, since this usually affects che 
viscosity exponentially, through an activation energy. However reducing 
the temperature from 19°c to —41°c only increased the modulus by 10% 
of the original reduction. Assuming that the coefficient of viscosity 
increased by a similar amount, this indicates an activation energy of only 
about 200 cals/mole, which is far below that normally encountered in 
viscous flow in solids, for which the activation energy is usually in the 
range 20 000-40 000 cals/mole. It therefore seems most unlikely that 


Zener’s explanation is correct. 


i00 HOURS 
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Both the other theories involve dislocations. Eshelby (1947) suggested 
that dislocations lie in potential energy ‘troughs’, and that these dis- 
locations oscillate (elastically) under the application of small alternating 
stresses, thus increasing the macroscopic strain and apparently reducing the 
modulus. However Eshelby attributed the whole of the reduction to the 
movement of edge dislocations, while in the opinion of Mott (1952) it is the 
motion of the relatively more mobile screw components of dislocation loops 
which is responsible. It is worth noting that the evidence indicates that 
it is probably an elastic oscillation of the dislocations which is occurring, 
rather than any form of hysteresis mechanism whereby dislocations * snap 
loose ’ at some small stress, such as that suggested by Nowick (1951) in his 
work on internal friction. This follows since an increase in strain amplitude 
by a factor of four did not affect the modulus, nor was the recovery speeded 
up by continuous vibration of the specimen. 

The Recovery Process : If Mott’s theory is correct, then on this model the 
recovery caused by low temperature annealing is consistent with his 
explanation (1951) of the cross-slip observed by Cahn (1951) in aluminium. 
Mott suggests that this is a consequence of the behaviour of two dislocation 
loops generated on different planes. When the screw component of one 
loop is above and parallel to the screw component of opposite sign of the 
other, thermal energy may allow the two loops to join up into a single loop, 
half in one plane and half in the other, joined together by two short edge 
dislocations. Another possible way in which annealing could cause 
recovery is by promoting the diffusion of solute atoms to ‘lock’ the 
dislocations, in the manner explained by Cottrell (1948). In either case 
there will be an activation energy associated with the recovery process, 
which will either be a measure of the energy barrier separating the two 
screw dislocations on Mott’s model, or the activation energy for the 
diffusion of the solute in copper. For the moment let us assume that it is 
the former mechanism which is operative here, and examine the experi- 
mental results on this basis. In order to obtain a value for the activation 
energy it is necessary to make some assumptions, which are : 


1. Only a given number, N, of screw dislocations lie in relatively 
favourable positions from which they can be moved by thermal activation. 

2. The rate of movement of dislocations with a given activation energy 
is proportional to the number still to be moved. 

3. The recovery in # at any time is proportional to the number of 
dislocations which have moved. 


4. The activation energies are distributed in a Gaussian manner about a 
mean. 


On the basis of these four assumptions equations are obtained in the 
Appendix from which both the mean activation energy, H, and the 
standard deviation, o, of the distribution may be determined from the 
experimental results. It is there shown that if the modulus after working 
is #,, after time ¢ is H and after the final prolonged low temperature 
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anneal is #,, then a good approximation to the relation between E and tis 


BE, .(-t  /-F 
Ba 01s | er |} 


where T is the absolute temperature, R the gas constant and 7, a constant. 
A simple way of obtaining H from the experimental data is to measure the 
timet, taken to reach a given degree of recovery at an annealing temperature 
T,, and to plot log, t, against 1/7',. Several such eraphs may be obtained 
for different degrees of recovery, which should be a number of parallel 
straight lines whose slope is proportional to H. | 

If the diffusion mechanism is the cause of recovery we should expect 
that the rate of diffusion of one solute element would dominate those of 
any of the others, so that all the dislocations would be anchored by this one 


Fig. 4 
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Time-temperature relation for a given amount of recovery : 1% extension. 


_O.F.H.C. copper. 
element. Hence a similar plot would again yield a series of parallel straight 
lines with a gradient proportional to the activation energy for diffusion of 
the solute in copper. 

One of the main differences between the two theories is that one would 
expect that the distribution of activation energy would be fairly wide if it is 
the thermal activation of the dislocations themselves which governs 
recovery, while it would be fairly narrow if it is the diffusion mechanism 
which is operative. Thus if we can determine the standard deviation, it 
should help to distinguish the two. 
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Interpretation of Results: Figures 4-6 show the graphs of logyot, 
plotted against 1/7’, for different recovery levels for the O.F.H.C. copper 
bars extended 1% and 8-4%, and for the spec. pure bars extended 1%. 
Examination of these figures shows that the points are scattered about 
straight lines, although the linearity in the case of fig. 6 is doubtful, since 
there are only three points on each line. The lines which have been drawn 
to fit the points have been calculated statistically as the best regression 
lines, assuming that all the error lies in log, t; and none in 1/7’, although 
had the opposite assumption been made the differences in gradient would 
have been small. 

Fig. 5 
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Time-temperature relation for a given amount of recovery : 8-4°% extension. 
O.F.H.C. copper. 


The standard deviation of the slope of each line was found by normal 
statistical methods and the 95% confidence zone found in the usual way. 
The mean values of these gradients (expressed as activation energies) 
together with these limits are given in table 1 for all the results shown in 
figs. 4-6, except that the 95° confidence zone for the spec. pure copper is 
so large, owing to the small number of results, that it is not worth quoting, 
and the 50° zone is given instead. 

It will be seen that there is a considerable difference between the 
activation energies for recovery for the O.F.H.C. copper extended by 1% 
and by 8:-4%. This suggests that the cause of recovery is more likely to be 
the annihilation of dislocations themselves rather than the diffusion of 
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solute atoms, since it is unlikely that the diffusion activation energy would 
change so appreciably after different degrees of working. It will also be 
shown below that the standard deviation of the distribution is such as to 
confirm this view, so that it will henceforth be assumed that recovery is 
caused by the elimination of the mobile screw dislocations. 


Fig. 6 
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Time-temperature relation for a given amount of recovery : 1% extension. 
Spec. pure copper. 


From table 1 a closer estimate can be made of H, by taking the variation 
in activation energy for the different recovery levels to be random, and 
weighting each figure according to its reliability. It should be realized 
that we have not obtained a number of independent estimates of H, which 


Table 1. Activation Energies for Recovery 


1% Extension 8-4% Extension 1% Extension 
E,—E ORE CS Cu O-F.H.C. -Cu Spec. Pure Cu 
Ei; = ———_. -of>asaq_~_~—_— 
; H(cals/mole) | 95% zone | H(cals/mole) | 95% zone | H(cals/mole) | 50% zone 
Ts 24 500 +21% — — — — 
6% 25 300 +15% 21 500 +18:5% = = 
5% 25 500 +13% 21 200 +4-0% 19 900 +11% 
4% 24 800 +14-5% 22 000 1.6-4% 21 700 £13% 
3% 25 200 +7:7% 20 200 aE 249, 22 100 +12% 
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may be used to determine the probable limits within which H lies, since 
errors which affect one value tend to affect all. It is in fact difficult to 
obtain any precise limits to set on 7, but they should not differ very much 
from the closest limits given in table 1. We may therefore estimate the 
95°, confidence zones, in round figures to be, 


1% extension O.F.H.C. Copper : 25 000 cals/mole+ 10%, 
84°, extension O.F.H.C. Copper : 21 500 cals/mole+-5%, 


while all that may be said about the spec. pure copper is that there is a 
50 : 50 chance that Af lies within the zone 21 000 cals/mole+ 12%. 

It should be noticed that if these confidence limits are accepted, and we 
regard the final mean figures as though they had been derived 
experimentally as the slopes of lines through 7 points (as is the case with 
most of the lines of figs. 4 and 5), the chance that the difference between 
the two values of H is not statistically significant is less than 1 in 1000, 

As mentioned above, the standard deviation, o, of the distribution of 
activation energies may be estimated from the gradients of the experi- 
mental recovery curves, figs. 1-3, by a method given in the Appendix. It 
is there shown that the values of o are : 


1% extension O.F.H.C. Copper: 1800-300 cals/mole, 
8:-4°, extension O.F.H.C. Copper : 2200-200 cals/mole, 


the limits placed on the values being the 95°, probability zones. This 
fairly wide distribution confirms the view that it is not the diffusion of 
solute elements which causes recovery. 

Summing up, we conclude that the most satisfactory explanation of the 
drop in modulus and subsequent recovery is that suggested by Mott (1952) 
viz. the application of an alternating stress causes relatively mobile 
screw dislocations to oscillate in their potential energy troughs, giving 
an extra contribution to the macroscopic strain which is entirely elastic ; 
low temperature annealing provides sufficient thermal energy to enable 
suitably placed screw dislocations lying in different planes to annihilate 
one another in pairs. From the results obtained it appears that the mean 
activation energy for recovery decreases as the cold-work increases, but 
that the standard deviation of the energy distribution increases. There 
is also some evidence that this activation energy is less for the recovery of 
spectrographically pure copper than for O.F.H.C. copper, for the same 
amount of working, but this is not certain. 
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PAE ee) INCL LX 


THE RELATION BETWEEN RECOVERY AND ANNEALING TIME 


The assumptions on which this relation will be derived have been given 
above in § 4. Now suppose that immediately after working, the number 
of dislocations contributing to the reduction in modulus with energies 
between H and H+dH isdn. Thenthe number of these dislocations which 
have moved after time ¢ is 6 (l—exp (—+#/r)), where it is assumed that 
T=7, exp (H/RT). Since the recovery in £ is proportional to the total 
number of dislocations which have moved, irrespective of their energies, 
the relation between # and ¢ is 


Soh ee re t H d i 
Petes We —exp| —“exp ~ Pr han, we ML) 


where JN is the total number of movable dislocations. If the activation 
energies are distributed according to the Gaussian law, the relation between 


nN and jel 18 n 7 i [ies i (HH) rT 
Nov (20) JF P GA , 


where H and o are the mean and the standard deviation, respectively. 
For any given values of H and o a good approximation to the solution of 
eqn. (1) may be made using a step-by-step calculation. The activation 
energy distribution is divided into equal zones, 1000 cals/mole wide, and 
the proportion of N having activation energies in each zone is found from 
tables. Then, taking H as constant in each zone, eqn. (1) is expressed as 
a summation of terms, each of which may be computed in terms of ¢, 7’ 
and the parameter 7). This computation has been made for ten different 
conditions, viz. a mean activation energy of 25 000 cals/mole with five 
standard deviations (sc=0, 500, 1000, 2000, 3000, cals/mole), for each of 
two annealing temperatures (7’=323 and 373°k). In order to test whether 
a zone 1000 cals/mole wide was sufficiently narrow one computation was 
made using zones only 500 cals/mole wide (c=1000 cals/mole and 
7T—373°K) ; no significant difference from the corresponding computation 
using wider zones was found. The results of the ten computations are 
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shown graphically in figs. 7 and 8 in which (E,—E)|(H,—E,) is plotted 


against (log, 9!—log ;979)- ' a 
The Aitene method shows how to obtain recovery curves given the 
distribution of activation energies: the practical problem is to find the 


distribution given the recovery curves. 


1. Determination of H. 


As a first approximation take H as constant in eqn.. (1) and equal to H, 


then E-E, Lf exp (exp (— aa) 
Tigo nti \ orm RT 
Now consider two recovery curves at different annealing temperatures 7, 


and 7',. Let the times taken to reach the same arbitrary level of recovery 
be ¢, and t, in the two cases. Then 


fo k HS (eae ea ) 
SRe re aT 
Applying this result to the theoretical recovery curves shown in figs. 7 and 


8 at different recovery levels and different standard deviations we obtain 
the values of H shown in table 1. 


Table 1. Values of H from Theoretical Recovery Curves 


H (cals/mole) 


o=2000 cals/mole o=3000 cals/mole 


24 000 23 400 
24 500 24 200 
25 000 25 000 
25 500 25 700 
26 LOO 26 400 


Now in this particular case the true value of is 25 000 cals/mole, so 
that table 1 shows that the first approximation taken above is accurate 
enough to give a reasonable figure for H, particularly in the region of 50°% 
recovery. 


2. Determination of o 

The value of the standard deviation of the distribution may be obtained 
from measurements of the gradients of the recovery curves at their points 
of inflexion. It will be seen from figs. 7 and 8 that this gradient is fairly 
sensitive to the value of oc, while some further computations showed that 
it is insensitive to the value chosen for H. The gradients of the experi- 
mental curves are shown in table 2, allowance being made for the fact that 
the first measured value of the modulus after working was not the true 
value of H,. 
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It is clear from this table that the experimental error outweighs any 
differences due to different annealing temperatures, so that the results are 
best treated as being random variations about a mean, corresponding to a 
mean annealing temperature of about 75°c. Comparing this value with 


Table 2. Gradients of Experimental Curves at Points of Inflexion 


1% Extension 8:4% Extension 
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Theoretical recovery curves assuming Gaussian distribution of activation energies. 
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the mean of the two gradients of the theoretical curves for different values 
of o, which mean is assumed to correspond roughly to an annealing 
temperature of 75°c, by interpolation the standard deviations of the 
distribution were found to be : 
1% extension 1800-+300 cals/mole, 
8-4% extension 2200+ 200 cals/mole. 


Thus both the parameters of the distribution may be obtained, and it has 
been shown that one depends but little on the other. 


3. Value of 7 
Comparison of the experimental curves with the corresponding 
theoretical curves, at 50% recovery, gives the values of 7) as approxi- 
mately 5x 10~™ seconds for the 1°% extension and 10~1! seconds for the 
8:-4°% extension. It is of interest that the period of one atomic vibration 
for copper is about 1-5 10~18 seconds, since this will be the shortest time 
for a dislocation to move one atomic distance. 
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ABSTRACT 


The adsorption of argon, nitrogen, and oxygen on an aluminium foil 
has been determined by a gravimetric method. The purpose of the 
investigation was to obtain data at high saturations for comparison 
with the Frenkel—-Halsey—Hill equation describing multilayer adsorption. 
The isotherms obtained indicate that this equation is obeyed closely 
for the three gases. Apparent deviations in the case of argon and oxygen 
can be ascribed to anomalous packing in the first few atomic layers and 
allowance can be made for this effect. The verification is discussed in 
terms of the theoretical basis for the equation. Data is also given for 
the cross sections of the argon and oxygen molecules in the first atomic 
layer. It is suggested that these cross sections are considerably greater 
than values calculated from the corresponding bulk density. The magni- 
tude of the standard nitrogen cross section is also discussed. 


$1. INTRODUCTION 

From a theoretical point of view, there are two regions of an adsorption 
isotherm which are of particular interest. These regions comprise, firstly, 
that at very low saturation, where the first few atomic layers are being 
formed, and secondly, the region at very high saturation, where the 
transition from multi-molecular layers to the liquid phase is taking 
place. Although there exists no comprehensive theory of adsorption 
describing the growth of the adsorbed layers over the entire pressure 
range, important theoretical calculations have been made which appear 
to be satisfactory within the two regions mentioned above. The physical 
difference between these two regions is best described in terms of the 
binding forces involved. At low saturations, the forces between adsorbed 
molecule and the adsorbing wall dominate the process, while near 
saturation, the forces between the molecules in the adsorbed layers 
become the most important. 

The formation of the first few adsorbed layers has been described by 
a number of theories, of which the most successful is that of Brunauer, 
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Emmet and Teller (1938); however, it is generally admitted that a 
number of crude approximations must be made in order to solve the 
equations describing the adsorption process. For example, it is assumed 
that the number of molecules adsorbed over a given site is independent 
of the numbers on adjacent sites. In other words, lateral interactions 
between molecules in the adsorbed layers are neglected. Furthermore, 
the energy of adsorption of all but the first layer is assumed to be equal 
to the energy of liquefaction, so that only the first layer is considered to 
be significantly different from the liquid. 

A number of attempts have been made to introduce refinements into 
the B.E.T. theory (e.g. Brunauer, Deming, Deming and Teller 1940, 
Cassel 1944), but a satisfactory solution for the complete adsorption 
isotherm has not been obtained because of the mathematical difficulties 
encountered when a more realistic model is assumed (Hill 1952). 

As a result of these difficulties, the region of high saturation is of 
special interest. Here the number of layers is high enough for anomalous 
properties of the first layer to be relatively unimportant. It is therefore 
possible to make predictions of the form of the isotherm without having 
recourse to the type of approximation necessary for any solution of the 
low saturation problem. Calculations for this region have been made 
by Frenkel (1946), Halsey (1948), Hill (1949), and McMillan and Teller 
(1951). They suggested that, approaching saturation, the isotherm 
should follow the relation . 


K 
—np/py= =» ea egret rh wy «ah 
where v is the number of layers adsorbed at a pressure p, and py is the 
saturated vapour pressure. K is a constant of proportionality determined 
by molecular forces, as will be discussed later. Actually, Halsey proposed 
a more general form 


K 
in P| Pyo= pe ’ 


without specifying the index s explicitly. 

Equation (1) has a simple physical basis. It is derived by assuming 
that the van der Waals forces between an adsorbed atom and the 
adsorbing wall fall off as the cube of the distance. The specific form 
of (1) results from equating the difference in the differential Gibbs free 
energy of the adsorbed phase and the bulk liquid, p»—y,=RT In P| Pos 
to the corresponding change in potential energy due to the van der Waals 
forces. It is possible to allow for the influence of the interactions within 
the adsorbed layer as well as the attraction due to the wall. The 
calculation of Hill (see McMillan and Teller 1951) gives the following 
form for the constant of proportionality in (1) : 
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where p is the molecular number density, d,;, the equilibrium distance 
between a pair of molecules of types i and j, €4 18 the potential energy 
of such a pair at the distance d,;, and o is the effective molecular cross 
section. The subscripts 1 and 2 refer to the adsorbate and adsorbent 
respectively. McMillan and Teller have shown that this constant is 


approximately 
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It is therefore possible to correlate the magnitude of adsorption at 
high saturations with the molecular forces involved. 

Further approximation enables us to express K in terms of the measured 
energies of adsorption. McMillan and Teller have shown that ¢,,—¢,, is 
of the order of (#,—£,)/N, where H, is the energy of adsorption of the 
first layer, H, is the energy of liquefaction, and N is Avagadro’s number. 
(H,—£)) is a quantity appearing explicitly in the equation of Brunauer, 
Emmet, and Teller and can be derived from the intercept of the B.E.T. 
plot of p/(~)>—p)V against p/p». (V is the volume adsorbed at the pressure 
p.) Thus we take for our final approximation : 


7 


= ort | 


Owing to the fundamental basis and essential simplicity of this relation, 
experimental tests of its validity are of particular interest. The data 
usually quoted as supporting these conclusions are those of Harkins and 
Jura (1944) who measured the adsorption of nitrogen and of water on anatase 
(powdered TiO,). These data showed that eqns. (1) and (3) are followed 
quite closely, although Halsey (1948) deduces that an exponent of v 
equal to 2-67 instead of 3 describes the data more accurately. Halsey 
considers this deviation from the cubic law to be important and has 
asserted that the empirical value of the exponent of v in eqn. (1) “has no 
simple relation to the power in the expression for London forces’. From 
the work of Drain and Morrison (1952), it appears that the exponent may 
even be as low as 2 for argon adsorbed on rutile. 

In the case of more complicated systems, it is evident from the work of 
Boyd and Livingston (1942) that for adsorption of water on graphite, and 
propyl] alcohol on anatase and barium sulphate, large deviations occur from 
the Frenkel-Halsey—Hill relation. This is to be expected, since the deriva- 
tion of eqn. (1) is only valid for adsorbed systems with spherical 
fields and this simplification may not be satisfactory for complicated 
systems where orientation effects in the surface layer are possible. The 
data of Boyd and Livingston are insufficient for an exact determination of 
the deviation. 

There are two reasons why further work on this high saturation region 
was desirable. In the first place, it is not easy to determine the number of 
adsorbed layers at high saturations using the volumetric method of previous 
workers, since the dead space corrections become large and deviations 
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from the ideal gas law make them difficult to calculate. Secondly, it is 
desirable to observe adsorption on surfaces which are better defined than 
that of a fine powder such as anatase, in which capillary condensation may 
produce serious errors, especially near saturation. McMillan (1947) has 
shown that many of the deviations from the ideal B.E.T. isotherm can be 
explained quantitatively in terms of surface heterogeneity, thus emphasizing 
the importance of the character of the surface. 

For these reasons a determination of the adsorption at high saturations 
for a number of gases has been undertaken. The measurements were 
carried out using nitrogen, argon, and oxygen as adsorbate gases and 
aluminium foil of low roughness factor as an adsorbent. Measurements 
have also been made with helium and these results will be described in the 
next paper (Bowers 1953). Due to the unique properties of helium at low 
temperatures, its adsorption is of interest for reasons beyond those being 
considered in this paper. 

The adsorption was determined gravimetrically, using a microbalance 
similar to that described by Gulbransen (1944) and more recently by 
Rhodin (1950). The principle of the method is simply that of determining 
the increase in weight of the specimen as adsorbed layers form on it. This 
method has the important advantage that at high saturations there are 
no corrections to be applied to the observed weight of the adsorbed layer. 
Furthermore, it is possible to use a relatively small specimen with a small 
roughness factor for this determination. 


§ 2, METHOD 


The microbalance used in this determination was very similar in design 
to that described by Rhodin (1950), having in particular the same sym- 
metrical construction in order to avoid buoyancy corrections. The 
apparatus was primarily designed for experiments on helium adsorption at 
liquid helium temperatures, but by filling the liquid helium dewars with 
liquid argon, nitrogen, or oxygen, the specimen and counterweight could 
readily be maintained at various temperatures between 77 and 90°K, 
determined by the boiling points of the respective gases. The arrangement 
is shown schematically in fig. 1. The specimen of aluminium foil and the 
aluminium counterweight, balanced to 10~>g, are supported from the arm 
of the balance by Kovar wires of 0-002 in. diameter. They are suspended 
inside tubes of } in. diameter which lead into two separate dewar vessels 
as shown in fig. I. 

The specimen of aluminium foil was wound into a loose coil: it had a 
geometrical area of 306 cm’ and weighed } gram. The counterweight was 
cylindrical and had a surface area of 2cm?. Thus the adsorption 
measured refers to the difference in area between specimen and counter- 
weight equal to 304cm*. The specimen was degassed initially at 350°c 
for 12 hours, the pressure in the specimen chamber being 2 10-® mm of 
mercury. This was done with the specimen in its final position in the 
cryostat and subsequent contamination was avoided by keeping the 
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chamber evacuated to 2x 10-6 mm between experiments. Presumably, 
there is a layer of chemisorbed oxygen on the aluminium, which is not 
removed in the degassing. 

Elaborate shock mounting was used in order to prevent vibration 
reaching the balance system. Much of the complexity in the apparatus is 
necessary only for the adsorption work with liquid helium ; accordingly, 


ADSORPTION 
“~~ PRESSURE p 


FOIL 


SOLID 
SPECIMEN COUNTERWEIGHT 


LIQUID 
HELIUM 
ARGON, NITROGEN, OXYGEN 


Schematic diagram of apparatus. 


the details of the microbalance and associated cryogenic equipment will not 
be described here. There is one aspect of the arrangement which differs 
from similar apparatus of other workers. When working with nitrogen, 
argon, and oxygen, it was found desirable to place these liquids not only in 
the dewars represented in fig. 1, but also in dewars which enclosed the ones 
shown. These outer dewars are essential for providing nitrogen temper- 
ature shielding of the inner dewars when liquid helium is being used, but 
their use with the other gases results in a quiescent liquid in the inner 
vessels and this reduces spurious movement of the balance, . The pumping 
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tube shown in fig. 1 is only used with helium, since delayed boiling in the 
other liquefied gases renders precise temperature control impossible 
during pumping. 
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Nitrogen isotherm at 77-4°K. 


§3. REsuLts 
Adsorption of Nitrogen 
In fig. 2 is shown an isotherm for nitrogen at 77-4°K and in fig. 3 the 
corresponding B.E.T. plot (Brunauer, Emmet, and Teller 1938) which 
indicates a satisfactory agreement between theory and experiment between 
5 and 30% saturation. The effective area of the foil determined from this 
B.E.T. diagram was 388 cm? as compared with the geometrical area of 
304cm?. This implies a surface roughness factor of 1:3. In the calcu- 
lation of the surface area, the cross section of the nitrogen molecule was 
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taken to be 16-24? (Emmet 1948, Harkins 1952, p. 227). Thus the 
specimen of aluminium foil could be considered relatively smooth from the 
point of view of adsorption, since this roughness factor is small. From the 
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B.E.T. plot it was calculated that a completed monolayer on the specimen 
weighed 11:1x10-$g; the number of layers could therefore be deter- 
mined to a precision of at least one-tenth of a layer. 
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We shall now turn to an analysis of these data. In fig. 4 we show a 
eraph of loglog p/p») against log (mass adsorbed). This graph should be a 
straight line, the slope giving the exponent of v in the Frenkel-Halsey—Hill 
relation, assuming the number of layers is proportional to the mass 
adsorbed. The slope of this straight line was found to be 3-0 as predicted 
by Frenkel and Hill. One notices that there is one point, corresponding 
to 99-3% saturation, which lies above the line, suggesting that in this 
very high saturation region the number of layers is increasing at a more 
rapid rate than given by the cubic relation. However, since the isotherm 
is extremely steep in this region, one cannot say with certainty that this 
deviation does not arise from experimental error. Nevertheless, it does 
seem possible that at the highest saturations the number of layers is 
indeed increasing faster than the cubic law would require. It is difficult 
to specify the exact meaning of such an increase, because even though the 
roughness factor of the aluminium foil is small, some form of capillary 
condensation may be influencing the total weight of the specimen so 
near the saturation point. It is probably significant that one can see this 
tendency to a more rapid growth in the number of layers at a lower satura- 
tion, of the order of 96°%, in the data of Harkins and Jura on anatase, 
where capillary condensation is more likely to occur. Apart from capillary 
effects, the influence of surface tension and gravity on the adsorbed film 
may have to be reconsidered at these high coverages. 

In fig. 5 we have plotted log p/p, against 1/v? in order to determine K, 
the constant of proportionality in eqn. (1). The agreement between 
experiment and theory is most striking when plotted in this manner. The 
value of the constant K for nitrogen was found to be 3-10. 

We may compare the empirical value with the approximate eqn. (3a). 
Substituting #,—H,;—800 cal/mole (Brunauer 1945, p. 157) in this 
equation, we obtain K~2-7 at 77°K. This is very close to the observed 
value of 3-1 mentioned above. This degree of agreement may be fortui- 
tous, since an order of magnitude is all that McMillan and Teller claim for 
their method of estimating K. 

Summing up the results for nitrogen, it seems that the Frenkel—Halsey— 
Hill relation is followed very closely; the proportionality between 
log p/p) and 1/r3 is established and the constant K has the correct order 
of magnitude. 


The Adsorption of Oxygen 
The adsorption of oxygen has been measured at 77-4°K and 89-7°K. 
The two isotherms were found to be similar in character, and in fig. 6 we 
show the data for 77-4°x. Adsorption and desorption points are shown in 
this figure. There is no evidence of significant hysteresis in the data. 
This is to be expected, since hysteresis in adsorption is usnally ascribed to 


cracks and pores in the surface of the adsorbent (Zsigmondy 1911, McBain 
1935). 
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The B.E.T. plot was found to be a straight line between 5% and 25% 
saturation. The weight of the first layer was found to be 10-lpe@ at 
77-4K and 9-05 wg at 89-7°x. Assuming the area for adsorption is 
388 cm*, as determined from the nitrogen isotherm, these values imply 
a cross section for the oxygen molecule in the first layer which is greater 
than that calculated from the bulk liquid density (Brunauer 1945, p. 287). 
The area per molecule in the first layer was found to be 20-3 a2 at 77-4°K 
and 22-7 A? at 89-7°K as compared with a value of 14-2 4? calculated from 
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the bulk liquid density at 89°x. The existence of anomalous densities in 
the first few layers is well known in the adsorption of gases (Beebe, 
Beckwith, and Honig 1945, Davis, DeWitt, and Emmet 1947, Emmet 
1948), although the values for the cross sections found in this experiment 
for oxygen seem higher than those of other workers. 

The existence of the anomalous first layer introduces a further problem 
in the comparison of our results with the Frenkel-Halsey—Hill equation. 
We must allow for the varying density of the layers in computing v, the 
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number of layers, from our measured values of the mass adsorbed. In the 
case of nitrogen, previous work has shown that the first layer is not 
anomalous, so that it has become a standard gas for area determinations 
(Harkins and Jura 1944); this means that we could take the number of 
layers as proportional to the mass adsorbed. It is instructive to see how 
the variation in layer density affects the evaluation of our data, because 
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this factor may be responsible for app: iati : 

ee ees p for apparent deviations from the Frenkel-— 
In fig. 7 we show a graph of log v against loglog p/P for data obtained 

at 89:7°K. In computing the points for curve I, we have assumed that all 

the layers have the same density as the first layer. The inverse slope of 

this line (which should equal the exponent of v in eqn. (1)) was found to be 

2-78. The points on curve II have been computed, assuming that only 
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the first layer has the anomalous density determined from the B.E.T. 
plot, while all the rest have the bulk liquid density. The inverse slope of 
this line is 3-34. It is quite clear that the assumptions underlying the 
computation of these two curves are unrealistic ; neither can be accepted 
physically. However, the errors involved are in opposite directions and | 
the real physical picture must lie between the two assumptions. The 
situation in the adsorbed film must involve a gradual change from the low 
density of the first layer to the liquid density in the outermost layers. It 
is interesting to see that one assumption yields an exponent slightly greater 
than 3, while the second yields one which is slightly less than 3. Since a 
realistic distribution of density would give a line in between the two shown 
in fig. 7, an exponent of 3 cannot be far wrong. 


Fig. 7 
p 
98 95 90 85.807570 po 


OXYGEN AT 89.7 °K 


o ood 


NUMBER OF LAYERS 
a 


-l0dio Mo 


The value of K computed from the experimental results is 4-5 for points 
on. curve I and 2-7 for those on curve IJ. These values have the correct 
order of magnitude as calculated from eqn. (3a). At the highest saturations, 
greater than about 98%, the number of layers again appear to be greater 
than expected from the cubic law. 

The isotherm obtained at 77-4°K for oxygen presents features which are 
very similar to the 89-7°K isotherm. Plotting the results as for curves I 
and II of fig. 7, we obtain inverse slopes of 2:7 and 3-1 respectively. The 
corresponding values for K are 6-4 and 2-9. (In comparing these values 
for K, it should be remembered that K determines the cube of the number 
of layers ; hence the evaluation of K is subject to 3 times the errors in the 
basic isotherm.) 

Thus the results for oxygen exhibit two features. They indicate the 
existence of a packing in the first adsorbed layer which is significantly less 
than that calculated from the bulk density. They also indicate that if 
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this anomalous packing is taken into account in a realistic manner, the 
isotherm at high saturation is described by eqn. (1). 


The Adsorption of Argon 
Isotherms for argon have been determined at 77-4°K and 87-5°K. At 
the saturation vapour pressure, argon is solid at 77-4°K and liquid at 
87-5°K. The isotherm at 77-4°K is therefore the only one involving solid 
layers which we have investigated. 
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The data for 87-5°K are shown in fig. 8. This isotherm exhibits features 
which are similar to the oxygen isotherms. The weight of the first layers 
was calculated to be 10:3 ug. This implies a cross section at the argon 
atom in the first layer equal to 24-8 A?, as compared with a value of 14-3 A2 
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calculated from the bulk liquid density. This disparity is even greater 
than that observed for oxygen ; it is similar in magnitude to that abeoreed 
by Beebe, Beckwith, and Honig (1945) for the krypton monolayer. Thus 
in analyzing the data at high saturations we are faced with the same 
situation presented by the oxygen results because of the anomalous 
first layer. In fig. 9 we show a graph of log v plotted against loglog p/py, 
computing v in the same way as for fig. 7. The assumption of a uniform 
layer density equal to that in the first layer (curve I) gives a value of 2-78 
for the exponent of v in eqn. (1). The assumption of a layer density given 
by the bulk liquid density in all but the first layer (curve IT) gives 3-13 for 
this exponent. Again we may assume that the real physical situation is 
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one involving the gradual change throughout the first few layers from the 
anomalous first layer packing to that calculated from bulk liquid density. 
Such a physical model would lead to a curve between I and II of fig. 9 
and at the highest saturations would be approximately a straight line 
whose inverse slope must be close to 3-0. 

In fig. 10 we show the data plotted on the form of 1/v? against —log P| Po; 
v was calculated as for curve II in fig. 9. The purpose of this plot is merely 
to show how closely the highest saturation data is represented by the cubic 
law in spite of the fact that an exponent of 3-12 should provide the best 
general fit. It emphasizes the question of whether we should consider small 
deviations from the power 3 to be significant, as has been done by Halsey 


(1948). 
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The value of K was found to be 4-7 for the points on curve II of fig. 9 and 
2-2 for the points on curve I. These values are of the correct order as 
calculated from eqn. (3@). 

The isotherm for argon at 77-4°K is shown in fig. 11. The value taken for 
Py is that appropriate to the solid phase at 77-4°x. Considerable difficulty 
was encountered when measuring the adsorption at high saturations in this 
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case. It is not certain that equilibrium was finally attained ; heat 
transfer between the specimen and the gas seemed to be ver ait: 
The isotherm has two interesting aspects. Firstly, the cross sue of the 
argon atom in the first layer is again greater than one can calculate from 
the density of solid argon ; the figures are 17-9 A? and 12-9 A2. respectivel 

Secondly, the isotherm appears to meet the p/py=1 axis at 2 finite pale 
in the region of 3-5 layers. The adsorption at high saturations of ohn 
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layers is little understood, but the latter effect has been observed before 
with other gases, for example, carbon dioxide (Brunauer and Emmet 
1937) and krypton (Beebe, Beckwith, and Honig 1945). As the result of 
this effect, the data at the high saturations cannot be represented by a 
power law similar to eqn. (1). 

Frequently the isotherm for solid layers is plotted with a p, appropriate 
to the supercooled liquid (David, DeWitt, and Emmet 1947) which is 
larger than the solid value for a given temperature. In the case of argon 
at 77°K, the difference between the two saturated values is 15% (Clark, et 
al. 1951). However, if the data are plotted with respect to this liquid pp, 
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the highest saturation obtained experimentally is 85%, so that the data 
are insufficient for a test of eqn. (1). 


§4. DISCUSSION — 

The results presented in this paper suggest that eqn. (1) is accurate in 
describing the adsorption of nitrogen, oxygen, and argon at the highest 
saturations and at temperatures where the layers can be considered to be 
‘liquid’. The apparent deviation in the case of oxygen and argon may be 
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ascribed to anomalous structure of the first few layers. The only serious 
deviations appear in the region of 99% saturation ; however, it is difficult 
to determine the isotherm precisely so close to saturation. 

The relation verified is based particularly on the following assumptions : 

(a) The van der Waals forces involved vary as 1/r3 (London 1930) ; 

(b) The distance of an adsorbed molecule in. the outer layer from the wall 
may be taken proportional to the number of layers. 

For the interaction of a neutral molecule with the metallic wall, calcu- 
lations of Lennard-Jones (1932), Bardeen (1940), and Margenau and 
Pollard (1941) all yield a 1/7 dependence for the forces involved. They 
differ only in calculations of the size of the force. Until the A of eqn. (1) 
can be estimated more precisely than by eqn. (3a), we are not able to make 
a comparison between our data and the calculations mentioned. Prosen 
and Sachs (1942) have made a direct quantum mechanical calculation of 
the interaction between a neutral molecule and a metallic wall, using 
perturbation theory. They propose that the potential energy varies as 
log 7/r? when electron degeneracy is taken into account. Of course, the 
interaction between neutral molecules and the non-conducting adsorbed 
film still varies as 1/73. The introduction of this potential into an equation 
similar to (1) would lead to deviations from the cubic law if its influence were 
significant compared with interactions between adsorbate molecules. 
However, it is almost certain that the treatment of Prosen and Sachs is not 
valid for the distances from the wall found in multilayer adsorption and 
should only be applied to limited aspects of monolayer adsorption. It 
breaks down at short distances where the wave functions of neutral atom 
and metallic electrons overlap and also at large distances where the 
electron . coulomb interactions cannot be properly treated. Thus 
assumption (a) seems justified. 

Halsey (1948) has criticized assumption (b). He points out that one 
should not take the distance of an adsorbed molecule in the outer layer as 
proportional to the number of layers, because the actual distance involved 
is not continuously varying, but is a step function with values which are 
appropriate to integral numbers of layers. However, this view is based 
on a rather idealized picture of the adsorbed layers and the validity of the 
criticism may be less when we consider the statistical nature of what we call 
‘the number of layers’. The verification of eqn. (1) for the simple 
systems described in this paper indicates that assumption (b) does not 
introduce serious error. 

In the course of this work it has been necessary to determine the density 
of the first adsorbed layers for nitrogen, oxygen, and argon on the specimen 
of aluminium. Unfortunately, the values we give for the cross sections of 
argon and oxygen are relative to the nitrogen cross section. If the 
nitrogen first layer is normal, as appears to be the case from other workers 
(Harkins and Jura 1944), then the first layers of oxygen and argon are 
significantly less dense than the outer layers calculated from the 
appropriate bulk density. The differences in density we obtain. are larger 
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than the comparable values of other workers (Emmet 1948). Considering 
the work of Beebe, Beckwith, and Honig (1945) on krypton, which yielded 
a disparity of similar magnitude, it is clear that the first layer packing of 
the isotherm is not adequately understood. 

With regard to the nitrogen first layer, Davis, DeWitt, and Emmet 
(1947) state the following: -“...it is well to keep in mind the present 
status of area determination by nitrogen and to ask particularly whether 
or not it is possible that the cross sectional area assigned to the nitrogen 
molecule may be too low. The most convincing answer to this question 
will be obtained when and if someone measures the surface area of a metal 
foil or ribbon using nitrogen as an adsorbate. Only if one then obtains 
too low a roughness factor for such surfaces, will one be able to conclude 
with certainty that the usual value of 16-2 A? per molecule is too small. The 
experiment of Harkins and Jura (1944) is not absolutely conclusive.” 
The value of 1-3 for the roughness factor of the aluminium foil deter- 
mined in the present investigation seems reasonable and lends support 
to the assumption of the normal 16-2 A? for nitrogen. 
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ABSTRACT 


The adsorption of helium on an aluminium foil has been measured 
between 1-8°K and 4:2°K. The principal interest was in the form of an 
isotherm close to saturation. The maximum thickness observed for the 
adsorbed layer at the highest saturations (greater than 99:9°%) was 
twenty atomic layers below the A-point and ten layers above the )-point. 
These values for the unsaturated film are compared with existing data for 
the thickness of the saturated helium IT film. The data are also compared 
with the Frenkel—Halsey—Hill equation for multilayer adsorption and the 
calculations of Band on the adsorption of a Bose-Einstein gas. 


§1. IyTRODUCTION 


Ir is well known that when any solid surface is in contact with liquid 
helium II or its vapour, a thick film is formed which possesses unique 
transport properties (Daunt and Mendelssohn 1939). Estimates of the 
thickness of the film above the liquid vary from 50 to 150 atoms. Recent- 
ly attention has been drawn to the fact that even in the comparatively 
thin adsorbed layers, characteristic of pressures below the saturated 
vapour pressure, superfluidity is observed when the temperature is 
reduced below 2:19°K. Long and Meyer (1950) showed that this peculiar 
mobility could be observed whenever the number of layers exceeded the 
first statistical layer. More recently, considerable work has been done 
on the magnitude of the mobility to be found in these thin layers (Bowers, 
Brewer and Mendelssohn 1951, Long and Meyer 1952a,b). In these 
experiments one measures the mass transfer directly or its contribution 
to the heat flow through the film. These total flows involve a product of 
the linear velocity of the flow and the thickness of the film, as well as 
p,/P, (the relative concentration of normal and superfluid components 
in the flowing film). One of the striking features of the data available 
is that the total superflow, even at 95% saturation, is still very much 
smaller (less than 10°) than that in the saturated film. It is of consider- 
able interest to determine whether this is merely a consequence of the 
much smaller film thickness involved in the unsaturated case, or whether 
the other factors determining the flow are also different in the two cases. 
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In order to separate the variation of the velocity and p,,/ps from the 
measured transport property, it is necessary to know the variation of 
the thickness of this film with pressure. Essentially, this requires an 
accurate determination of the adsorption isotherm. There is available 
a considerable amount of data on the isotherm up to about 90% saturation 
(e.g., Keesom and Schweers 1941, Frederikse and Gorter 1950, Schaeffer, 
Smith and Wendell 1949, Strauss 1952) and, on the whole, the agreement 
between the results of different workers appears to be satisfactory. 
However, much less is known about the form of the isotherm at the 
highest saturation. The isotherms determined by Kistemaker (1947) 
and Long and Meyer (1949) are the only ones which give measured points 
close to saturation. Long and Meyer found that the thickness of the 
film is of the order of 7 to 8 layers at 90° saturation, as observed by the 
earlier workers. However, with helium II, the adsorption at saturations 
higher than 90°, was found to be anomalously high, the isotherm rising 
very rapidly close to saturation to a measured value of 160 atomic layers. 
This effect was not observed above the dA-point. In Kistemaker’s 
original paper, he calculated about 30 layers for a similar limit, but a 
re-evaluation of his data (Frederikse 1950) indicates a film several 
times thicker than the 30 layers stated. The observation by Long and 
Meyer and by Kistemaker of a very thick film close to, but definitely 
below, the saturation pressure is of considerable importance, since it has 
direct bearing on the problem of whether the thick helium IT film is 
formed from the vapour, or whether contact with the liquid is necessary 
for its formation (Mendelssohn 1946). If we assume that they are 
observing the formation of the thick helium II film, normally character- 
istic of saturation conditions, there is a large disparity between their values 
and the most recent direct determinations of the thickness of the saturated 
film (Burge and Jackson 1951, Henshaw and Jackson 1951). The latter 
workers have obtained thicknesses which are less than a third of the 160 
atoms recorded by Long and Meyer. 

The form of the helium isotherm at high saturation is of interest from 
another point of view. One may ask whether the Frenkel—-Halsey—Hill 
relation, described in the preceding paper (Bowers 1953), is valid for helium 
adsorption. This law was actually derived on a classical basis, using the 
1/r? van der Waals force. If the cubic law is found to be valid, it would 
mean that quantum phenomena influence the magnitude of the forces 
involved, but not their dependence on the effective distance. On the 
other hand, quantum effects such as Bose—Kinstein degeneracy, might 
dominate the picture to such a degree that the simple cubic law is no longer 
valid. The influence of the gravitational field on the Bose-Einstein 
condensation is also a factor that may influence the adsorption consider- 
ably (Lamb pea Nordsieck 1941, Becker 1950, Liebfried 1950, Halpern 
1952 a, b). Calculations on the influence of degeneracy on helium 
adsorption. have been made by Temperley (1949) and Band (1951). 

This paper describes a determination of the helium isotherm by a 
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gravimetric method. The work has been undertaken for the following 
reasons : : 

(1) To investigate the high saturation region in greater detail in order 
to obtain data which are sufficient for a test of the cubic law. 

(2) To investigate the disparity between the values reported by 
Kistemaker and by Long and Meyer for the thickness of the film and the 
comparable values due to Jackson and co-workers. 

(3) To provide a basis for a gravimetric determination of the thickness 
of the helium ITI film under saturated conditions. 


Fig. 1 
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Helium isotherm at 1-:829°k. The symbol a denotes the arbitrary starting 
point mentioned in the text. 


§ 2. MerHop 


The apparatus and specimen used in this experiment were the same as 
used for the nitrogen, oxygen, and argon work described in the preceding 
paper. The microbalance and associated low temperature system will be 
described in detail elsewhere. 

The two helium baths (see fig. 1, preceding paper) were pumped through 
a common pumping tube, keeping the same temperature in both. Further- 
more, care was taken to keep the helium levels in these two dewars at the 
same height. The vapour pressure of the helium in the dewars is the py of 
the isotherm ; the value of p is the pressure in the tubes containing the 
specimen and counterweight. The specimen and counterweight space 
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could be evacuated to 5X 10-7 mm Hg when liquid helium was in the 
dewars. During a determination of the isotherm, the vapour pressure Po 
in the helium dewars was held constant to 0-1 mm of oil (octoil-S). The 
pressure inside the specimen chamber was changed by adding purified 
helium gas and, as in the earlier experiments, the deflection of the balance 
beam was observed as the adsorbed layers increase the weight of the 
specimen. Two radiation shields, with orifices of 1 mm, were mounted 
above the specimen. Both shields were in good thermal contact with the 
outer helium bath. . 

In spite of the elaborate shock mounting of the apparatus, a slight 
residual disturbance was noticeable when the helium was being pumped. 
This presented a difficulty only when a zero reading was being taken with 
the specimen space evacuated ; the presence of even a small amount of gas 
in the specimen space provided enough damping of the balance for complete 
stability. Accordingly, all the isotherms below 4:2°K are measured from 
an arbitrary starting-point at about 15% saturation, and in computing 
the total number of adsorbed layers, it is necessary to allow for the gas 
already adsorbed at this initial saturation. Such a correction can be made 
quite accurately, using the data of other workers, because at 15% satur- 
ation we are in the region of the completed monolayer where the isotherm 
has a small slope. Since no difficulty of a similar nature was encountered 
when using liquid nitrogen, oxygen, and argon at their boiling points, the 
disturbance presumably originated in the pumping system. A zero 
reading was obtained for helium at 4:2°K, 


§ 3. RESULTS 


Isotherms have been determined at the following temperatures : 
1-644, 1-829, 2-050, 2-132, 3-012, 4-198°K. In fig. 1 is shown an isotherm 
for helium IT determined at 1-829°K. The left-hand ordinate gives the 
mass of helium adsorbed, while the right-hand ordinate gives the equiva- 
lent number of statistical layers. For layers other than the first two, the 
weight of the layer is calculated assuming a normal liquid density, taking the 
effective area of the specimen as 388 cm? as determined by the nitrogen 
isotherm described in the preceding paper. We allow for the fact that the 
density in the first and second layers is significantly higher than successive 
layers by the following method : we have assumed that the density in the 
first layer is 3-6 times that for liquid packing (Schaeffer, Smith and 
Wendell 1949, Long and Meyer 1949, Gorter and Frederikse 1950, Strauss 
1952), while the second layer is arbitrarily taken to be half as dense as the 
first. It can be seen that the right-hand ordinate in fig. | has a scale 
which allows for these closely-packed layers. 

In fig. 2 we have a similar isotherm at 2-132°K, including in this curve 
some values for the coverage obtained by desorption from 99-25% 
saturation. It can be seen that there is no significant hysteresis, as was 


the case for oxygen with the same specimen (see preceding paper, Bowers 
1953). 
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As a final example of a measured isotherm, we show in fig. 3 the data for 
4:197°K. Two isotherms have been obtained above the \-point, at 3-02°K 
and 4-20°K. It is much more difficult to obtain accurate isotherms in the 
helium Irange. The bubbling of the liquid helium I and the simultaneous 
difficulty in controlling temperature make the balance less stable. 
Furthermore, the total number of layers found above the A-point is con- 
siderably less than the number observed below the A-point for the same 
saturation. The maximum number of layers observed above the A-point 
was ten at 99-91%, saturation, for the 4-2° isotherm.* 


Fig. 2 
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Helium isotherm at 2:132°xK. 


§4. Discussion 

The results presented above will be considered from two points of view. 
First we shall consider their relation to the formation of a very thick film 
at saturation and, secondly, we shall consider to what extent the adsorp- 
tion at high saturation is compatible with a law of the type proposed by 
Frenkel, Halsey, and Hill outlined in the introduction to the preceding 
paper. With regard to the formation of a thick film at saturation in 
helium II, it can be seen that there is no sign of a 160-layer film even at 
saturations as high as 99-92°%. However, the isotherm is very steep in this 
region and it is difficult to put a definite upper limit on the number of layers 
which one could obtain, but in no case were films thicker than 20 atoms 
observed. Furthermore, it may be said provisionally that we have not 
been able to obtain films thicker than 40 or 50 layers even with the bulk 
liquid helium II present, providing the specimen is in contact with the 
vapour alone as opposed to being in direct contact with the liquid (see 
panes eiaves ieee MMMM gL iied IPE IST nian 4) cS es 


* See recent paper Jackson, L. C., and Henshaw, D. G., 1953, Phil. Mag., 
44, 14. 
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Mendelssohn 1946). (Work on the film see with the bulk liquid is 

ill in progress and requires further clarification. 
pies eae seem a contradict those of Kistemaker and of Long and 
Meyer, which yielded a film about 160 atoms thick near saturation. Our 
results are in better agreement with those of Jackson and co-workers. 
One explanation of the disparity with the results of Long and Meyer may 
be found in capillary condensation, since they used tightly-packed rouge 
as an adsorbent. However, the fact that Long and Meyer did not find a 
similar thick film above the A-point is evidence against capillary 
condensation. 


Fig. 3 
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There is another possible explanation for the disparity involving the 
dependence of the helium II film thickness on the height above the bulk 
liquid. It is well known that the thickness of the film varies rapidly with 
the height above the bulk liquid level for heights less than a few milli- 
meters (Burge and Jackson 1951). Any bulk liquid introduced into a 
tightly-packed powder, such as the rouge, would be dispersed throughout 
the rouge by the large surface tension forces in the very small channels, and 
if we imagine a helium IT film being formed from this dispersed liquid, it is 
clear that the film on the grains will have a thickness characteristic of a 
‘height ’ of the order of the grain size. Since this effective height is 
extremely small, we should expect the thickness to be considerably 
greater than the value obtained in experiments such as those of Burge and 
Jackson, where a height of order of the centimetres is involved. Thus, 
even if a correction could be made for the presence of bulk liquid in an 
adsorption experiment at saturation, it seems reasonable to expect a film 
thickness for the fine powder which is greater than that obtained in other 
determinations. This is, of course, largely an academic point, since such a 
correction for the presence of bulk liquid is almost impossible to make. 
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None of the arguments presented above can be used in order to explain 
the disparity between our values and those of Kistemaker. He measured 
the adsorption of helium in a glass system designed so as to virtually 
eliminate the possibility of capillary condensation. At the present time 
no reasonable explanation of this disparity presents itself. The thickness 
of the film very near saturation should be extremely sensitive to the 
temperature of the specimen with respect to its surroundings and the 
disparity may indicate that the temperature of the specimen was higher in 
our case. This is mentioned because it is definitely easier to shield the 
specimen in the volumetric methods of Kistemaker and of Long and Meyer 
than in the present investigation. However, considering the radiation 
shields used in the present work, this mechanism for the disparity remains 
unlikely. 

We have already mentioned in the introduction to this paper, the fact 
that the total superflow observed in these adsorbed layers is very much 
smaller than that in the saturated film. Our isotherms and those of 
previous workers indicate that this is largely due to the smaller thickness 
of the adsorbed layers as compared with the film at saturation. However, 
there is a spread of about 50 to 150 atoms in the existing values 
for the thickness of the saturated film, so that it is not possible to 
analyse this aspect of the isotherms in any detail. If the saturated film 
is 50 atoms thick, the reduced thickness of the adsorbed film at 95%, 
saturation would result in a direct reduction of superflow to about 20% 
of its value in the saturated film ; a saturated film 100 atoms thick would 
make this figure 10%. 

Leaving the question of the relation to the bulk liquid film, we shall now 
consider the over-all magnitude of the adsorption observed. Our data are 
in general agreement with the previous work of Gorter and Frederikse 
(1950), Kistemaker (1947), Long and Meyer (1949) (except as mentioned 
above), and Strauss (1953). Generally speaking, helium forms a larger 
number of layers above about 70°, saturation than one finds with other 
gases (cf. preceding paper, Bowers 1953) when the temperature is below the 
A-point. 

In order to compare our data with the equation of Frenkel—Halsey—Hill 
(eqn. (1) of the preceding paper) we have plotted in fig. 4a graph of log p/p 
against 1/v? for three different temperatures. As in the previous 
notation, v is the number of layers. It is clear that the data illustrated 
can be described by the relation : 

Deke 

SA eee 

Po 
where K varies with temperature, and this seems to be true of all the 
isotherms obtained. The very striking way in which the data at 3-02°K fit 
on a straight line is to some extent accidental, since this degree of agree- 
ment exceeds the experimental precision for isotherms above the A-point. 
However, taking errors above the A-point into account, it seems as though 
the relation is consistent with the data between 1-8°k and 4:2°xK. The 
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variation of the constant of proportionality K with 1/7 is shown in fig. 5. It 
must be remembered that the value of K determines the cube of the 
number of layers, hence its determination is subject to three times the 
errors inherent in the original measurement of the isotherm. The variation 
of K in the helium II region from 60 at 1-64°K to 24 at 2-13°K involves only 
a difference of 35° in a number of layers. Hence the points in fig. 5 
should be taken to illustrate the definite rise in K in going through the 


Fig. 4 
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A-point, but they are not accurate enough to give the exact functional 
ae of K ontemperature. (In fig. 5 a straight line has been drawn 
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In the preceding paper, estimates of K have been made using the 

approximation : 
7 
K = opp (Hi—H), 

where H, is the energy of adsorption in the first layer and H, is the energy 
of liquefaction. A difficulty arises in the application of this formula to 
helium adsorption because H,—E, varies considerably with coverage 
throughout the first layer (Keesom and Schweers 1941). If we obtain 
E iL, from the B.E.T. plot, as in the preceding paper, we obtain a value 
appropriate to the completed first layer which has an anomalously high 
density. The B.E.T. plot of Schaeffer, Smith, and Wendell (1945) yields 


Fig. 5 


The variation of K with 1/7. 


about 25 cal/mole for H,—H;. However, the data of Keesom and Schweers 
(1941) show that #,—H,is at least 80 cal/mole at the lowest coverages and 
one could justify taking even higher values on the basis of the anomalous 
high density of the completed first layer. Taking #,—H,=25 cal/mole, 
we obtain K ~3, while an assumption of 80 cal/mole yields K = 10. It is 
important to emphasize that these values for K represent no more than a — 
guess. Firstly, the value taken for #,—H; is arbitrary. Secondly, on 
more fundamental grounds, we cannot expect a formula derived classically 
to explain helium II adsorption which is probably influenced considerably 
by the Bose-Einstein degeneracy. Apart from statistics, we have allowed 
for the anomalous first layer in computing the total number of layers by 
assigning a particular high density to this layer, but we have not considered 
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the possibility of a cooperative effect between the first and subsequent 
layers which would render the treatment of Hill and Frenkel inapplicable. 

The experimental data are best described by a K=7 in the helium I 
region ; this is not too far from the estimates made above. In the helium 
II region it is clear that the experimental K is considerably larger than the 
estimates. The data show that K is not proportional to 1/7’, but this is 
not surprising because of the A-transition. Summing up the observations, 
it is found that the Frenkel—Halsey—Hill equation describes the variation 
of the number of layers with pressure quite well for a given temperature ; 
the variation with temperature is complicated by the A-transition which 
invalidates the application of their formula. 


Fig. 6 
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The solid lines are the isotherms determined in the present investigation. The 
broken lines represent the calculations of Band with n=2, 2:5 and 3; 
the highest broken line refers to n=2 and the lowest to n—3. 


Recently Band (1951) has provided another method for the comparison 
of our results with theory. He has calculated the form of an isotherm for 
a degenerate Bose—Hinstein system, taking various powers in the force 
between. the adsorbate gas and the solid adsorbent. In fig. 6 we domipare 
his predictions with our data. In this figure are shown six isotherms and 
three calculated curves of Band. The various dotted lines correspond to 
n=2, 2-5, and 3, where n is the exponent in the attractive field obentiay 
proportional to I/r". It appears that the curve with n—2 fits Thoda 
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most successfully. Band’s formula implies that the isotherm is independ- 
ent of temperature for a fixed value of n and clearly this is an approximation 
in the case of helium. 

It is difficult at the present time to choose between the equation of 
Frenkel—Halsey—Hill and the equation of Band as representing the most 
reasonable picture of helium adsorption. The agreement with the data 
in the case of the former suggests that a binding potential energy propor- 
tional to the inverse cube of the distance is applicable to helium at low 
temperatures, but the magnitude of these forces may be very different 
in the helium IT region from that which one could calculate classically. 
On the other hand the description given by Band, although based on the 
important degeneracy, has one feature which is difficult to understand 
physically. As mentioned above, the best fit between his formula and the 
helium IT data is given when a binding potential between the adsorbate 
gas molecule and the combined wall and helium film varies as the inverse 
square of the distance involved. It is not easy to see the significance of 
such a law of force (see discussion in preceding paper, Bowers 1953). Such 
a formulation may be compared with the calculation of Bijl, de Boer, and 
Michels (1941) on the influence of zero-point energy in a thick film of helium. 
However, this calculation has been criticized by Mott (1949) because their 
potential term in 1/r?, which results from the zero-point energy, would not 
exist if lateral interactions in the film were taken into account. Temperley 
(1949) has also calculated quantum-mechanically the binding energy of a 
gas molecule with a wall and adsorbed layer. He considered the approxi- 
mation of a very thick adsorbed layer. This calculation yields an energy 
with terms in 1/r? and 1/r+... However, the term in 1/r? is very much smaller 
than the term in 1/74, so that this cannot be considered an interpretation 
of the Band formula. 

The influence of the zero point energy is obviously important in calcu- 
lating the magnitude of the adsorption of helium, but one can see from the 
large change in K in the region of the A-point that zero-point energy alone 
cannot explain the disparity between a classical treatment of the problem, 
such as that of Frenkel and Hill, and the observed isotherms. It would be 
interesting to see if one could calculate its influence on K ; no calculation 
exists at the present time on this point. 


§5. CONCLUSION 


The isotherms obtained with helium possess the following features : 
first, there appears to be no evidence with helium II for a film greater than 
about 20 atoms thick below the saturation pressure. Secondly, the over-all 
magnitude appears to be best described by a relation derived from a van 
der Waals force varying as the inverse cube of distance. The results 
can also be described reasonably successfully by Band’s formula for the 
adsorption of a degenerate Bose-Einstein system. In this case it is 
necessary to assume that the force of interaction between a gas molecule 
and the wall varies as the inverse square of the distance. 
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ABSTRACT 


The various mechanisms for the conversion of absorbed y-ray energy to 
heat are surveyed. It is shown that in dielectric crystals electronic 
excitation of the constituent molecular groups, followed by dissociation 
and recombination, and ionization followed by recombination, are to a 
large extent responsible for the energy dissipation. The partial inhibition, 
at low temperatures, of the final conversion to heat is discussed, with 
particular reference to the y-ray heating method in low temperature 
calorimetry. 


In research at temperatures below about 1°k, a major problem is the 
conversion of readily measured ‘ magnetic temperatures ’ to true thermo- 
dynamic temperatures (Cooke 1949). In one essential stage of this 
conversion, namely, the experimental determination of heat capacities 
on the ‘ magnetic temperature scale’, it is necessary to add heat to a 
sample at an ascertainable rate. One feasible way of accomplishing this 
is exposure of the sample to a beam of gamma-rays (Kurti and Simon 
1938). It appears to be tacitly assumed in the literature (Kurti and 
Simon 1938) that the energy so absorbed by the sample ultimately 
manifests itself entirely, or principally, in the form of heat. This note 
draws attention to the fact that such an assumption is unjustified. 

It is necessary first to trace the course of the absorbed energy through 
various stages of its dissipation. When a material is exposed to gamma- 
rays (or x-rays), those photons which are absorbed produce swiftly 
moving electrons (Compton recoil electrons, photoelectrons, electron 
pairs), virtually all of which have kinetic energy very great compared to 
the binding energies of valence electrons in the material. Thus, at this 
‘ first stage ’, the absorbed energy is still in a form which may be thought 
of as having high thermodynamic potential. The fast electrons interact 
with the medium and are gradually decelerated and ultimately captured 
(in the course of deceleration other (secondary, etc.) electrons may be 
ejected, these being similar in effect to primary electrons). The energy 
loss of all such electrons, as long as their energy exceeds the lowest 
electronic excitation potential of the medium, has the character that 
electronically excited and ionized atoms and molecules are predominantly 
produced. Some energy is transferred from the moving electrons to 


vibrational modes of the lattice, but this process is highly improbable 
SE eee 
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(Platzman 1952), except for the slow electrons for which it is the only 
possibility ; hence the fraction of originally absorbed energy so transferred 
must be small, certainly less than 15°% and perhaps not in excess of a 
few per cent. Thus the ‘ direct > production of heat is comparatively 
unimportant. Most of the heat follows later from the electronic excitation 
and ionization, i.e. the absorbed energy, in being degraded from electron 
kinetic energy to heat, passes, for the most part, through the intermediate 
form of atomic or molecular potential energy. At this ‘ second stage * of 
dissipation, then, most of the energy has retained an essentially high 
thermodynamic potential. 

It will now be argued that at the conclusion of a third and ‘ final stage ’, 
a substantial fraction of the energy originally absorbed will persist in 
forms other than heat and will, indeed, at extremely low temperatures, be 
stable with respect to transformation into heat. 

The ionized atoms and molecules are, in most solid substances at 
ordinary temperatures, quickly neutralized through various types of 
recombination. The mechanisms here are perhaps not understood in all 
detail, but one does know that the recombination energy is usually 
converted to heat more or less directly (often through the agency of 
internal conversion; see below). At extremely low temperatures, 
however, the recombination is likely to be inhibited and much of the 
recombination energy stored in the substances (as free energy), for the 
ejected electrons may be trapped at remote sites, the energy released in’ 
the act of trapping being dissipated to heat but the remainder of the 
recombination energy being preserved because of the smallness of the 
mobility of the trapped electrons. (The ionization energy is partially 
stored in some crystals even at room temperatures ; consider the case of 
colour centres.) 

The fate of excited molecules is next considered. (Only electronic 
excitation energy is treated ; some vibrational excitation almost always 
accompanies both electronic. excitation and ionization (Franck—Condon 
principle), but this relatively small quantity of energy is quickly dissipated 
to heat.) Ifthe excited system were an isolated atom it would radiate the 
excitation energy, but the case of an atom coupled to other atoms in a 
lattice, or of a polyatomic molecule, is entirely different. Here radiation 
is comparatively improbable and may in general be disregarded (although 
in exceptional cases it may occur with high yield; consider cathodo- 
luminescent or scintillant materials). Instead, the two most likely 
consequences of electronic excitation are dissociation (including pre- 
dissociation) and internal conversion (Teller 1937, Franck and Sponer 1948), 
The latter process, very common in polyatomic molecules, results in a 
direct conversion of the excitation energy to molecular vibrational energy 
(and thence quickly to heat). The former, however, results in molecular 
fragments (which will always have some kinetic energy, this being quickly 
dissipated to heat) having considerable chemical free energy. . Such 
fragments may, in rare cases be stable, or they may react chemically with 
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other molecules in the vicinity, or, most commonly in condensed phases, 
they may recombine. This recombination, because it may be followed by 
internal conversion, therefore expedites ultimate dissipation of the 
excitation energy to heat. It is important to appreciate, however, that 
such recombination almost invariably has some activation energy, and that 
this activation energy is such that at temperatures below 1°K the 
recombination is completely inhibited. Consequently, that portion of the 
absorbed energy which is utilized to dissociate molecules is for the most 
part lost as far as heating is concerned. In experiments at extremely low 
temperatures the working substance is a paramagnetic crystal (e.g. 
various alums) and in almost all cases it contains polyatomic groups such 
as NH,, SO,, and H,O which have great probabilities for dissociation if 
electronically excited. Indeed, ionization of such molecular groups is 
also commonly followed by dissociation. 

It may be mentioned that in these crystals virtually no energy is 
transferred directly from the moving electrons to the energy levels of the 
paramagnetic ions which must ultimately acquire the thermal energy at 
very low temperatures. Only after the absorbed energy is dissipated to 
molecular and lattice oscillation does energy transfer to the ions occur. 

Quantitative information bearing on these complex phenomena is at the 
present time insufficient for the calculation of the fraction of absorbed 
energy which is stored and thus unconverted to heat, but an estimate for 
very low temperatures ranges from about 20 to 70%. The result of 
irradiation by alpha-particles would be approximately the same, although 
there is a possibility of slightly greater efficiency of conversion because of 
high transient ‘temperatures’ in the vicinity of the alpha-particle 
path. The conversion efficiency could be measured approximately by 
repeating the Kurti-Simon experiment, and at the same time computing 
the energy absorption by means of an absolute determination of the 
intensity of the radioactive source and the various absorption coefficients. 
Tt would, of course, be expected to differ for different substances. (It would 
not, in general, be 100% even at normal temperatures, and for certain 
substances would apparently exceed 100%.) 

The possible consequences of these considerations in interpretations of 
experiments at extremely low temperatures will not be analysed here. It 
may be mentioned, however, that the method of Kirti and Simon (1938) 
is to a first approximation protected from error by the fortunate fact that 
the heat absorbed from the gamma-ray beam is not calculated from the 
intensity of the source, but, rather, is deduced indirectly, and in such a way 
that the only assumption regarding heating is that the heat added per 
unit time be independent of time, and of temperature. However, the 
possibility that this heating, in view of the mechanisms discussed apovss 
may not be strictly independent of temperature in the range from 0-01 
to 1°k, might merit examination. (There exists ample evidence of 
non-linearity with time, and dependence on previous treatment, in the case 
of alkali halide crystals irradiated at room temperatures (Kstermann, 
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Leivo and Stern 1949).) It should also be borne in mind that the sample 
is altered chemically by the irradiation, paramagnetic atoms and radicals 
being a common product; moreover, if repeated experiments are 
performed on a sample without permitting its temperature to rise greatly 
in intermediate stages, this decomposition is cumulative. 
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ABSTRACT 


The effects due to delayed heat evolution that might accompany y-ray 
absorption in crystals are discussed. It is concluded that, at least in the 
experimental conditions and with the substances employed so far, y-ray 
heating is a reliable method for use in calorimetry in the ‘ magnetic’ 
temperature region. 


THE foregoing paper by Professor Platzman (1953) on the mechanism of 
y-ray absorption in dielectric crystals and its bearing on low temperature 
calorimetry offers a good opportunity of stating briefly the essential 
features and some practical aspects of the y-ray heating method that has 
been used extensively in the Clarendon Laboratory, and, in some experi- 
ments, at the Laboratoire du Grand Eléctroaimant, in Bellevue, Paris 
(Kurti and Simon 1935, 1938 a, Kurti, Lainé and Simon 1937, Cooke 1949, 
Bleaney 1950). Such a survey is particularly desirable since some doubts 
have been expressed about the reliability of the method (de Klerk, 
Steenland and Gorter 1948, 1949). 

When calorimetric experiments below 1°K (in the temperature range 
obtainable by adiabatic demagnetization) were first begun it was realized 
that conventional heating methods were unsuitable below about 0-1°k 
because of the poor thermal contact and thermal conductivity at 
these temperatures. Heating by absorption of radioactive radiations 
appeared to be a suitable alternative (Simon 1935, Kurti and Simon 1935) 
because this would give uniform heat generation throughout the specimen 
provided that the attenuation of the radiation intensity across the speci- 
men was small enough. This condition could be easily satisfied for 
y-rays of, say, 1 Mev, with specimens of about 1 cm thickness and a 
suitable geometrical arrangement of the y-ray sources. For instance, if 
the sources have to be disposed rather close to the specimen so as to get 
the required heat input, variation of intensity across the specimen can be 
considerably reduced by rotating the sources round the cryostat (e.g. 
Hull 1947). 

Two possible shortcomings of y-ray heating have had to be considered. 
- One was the difficulty of obtaining an absolute value for the heat gener- 
ated: this was overcome by calibration in a temperature range in which 
the thermal properties of the specimen were known from other experi- 
ments. The other and more serious difficulty was the possibility of an 
de) * # OotAiniinicated by the Authors. 
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incomplete conversion of absorbed y-rays into thermal energy or, more 
precisely, of the existence of a finite relaxation time for this energy 
transfer to the system of electron spins. It was realized, of course, that 
with our method of determining the absolute value of the rate of heat 
input no error could be caused by this so long as the relaxation time was 
either very short or very long compared with the appropriate phase of the 
experiment. But for intermediate relaxation times the error might have 
been considerable. However, no indications of such relaxation times have 
been found in the numerous experiments with y-ray heating. 

Let us first consider the case of a relaxation time comparable with a 
heating period during a specific heat measurement. This would show 
itself by a delay in the temperature rise at the beginning of the period and 
a delay in the establishment of a steady drift rate after the heating period. 
No such time effects were usually observed and, on the few occasions when 
they did occur they were found to be due to inhomogeneous heating by 
y-rays (Hunt 1952). It is of course very important to arrange the y-ray- 
sources in such a way as to produce a uniform incident intensity over the 
whole specimen. Failing this, considerable temperature differences may 
be established in the specimen ; their equalization results in time effects. 
The absence of time-effects shows that, at least for the substances used by 
us, the relaxation time for the y-ray conversion is either less than half a 
minute or more than about five minutes. 

Next we have to consider the possibility of a relaxation time of between 
five minutes and the duration of a calorimetric run (the time taken to 
cover the entire temperature range, from about 0-01°K to 1°K), that is, 
of the order of a few hours. In this case, the effective rate of heat evolution 
would be a function of time and would, in particular, depend on previous 
history. One would, for example, get a different value for the specific 
heat at a given temperature according to whether one measured it on a 
‘ virgin ” specimen, i.e. not exposed to y-rays in the preceding hour or so, 
or whether this particular determination had been preceded by y-ray 
heating. Since no discrepancies of this type were found in our experi- 
ments relaxation times of this order may also be ruled out. 

There is also the possibility, as Professor Platzman points out, of a 
temperature-dependent relaxation time, very long in the lowest temper- 
ature range and rather short in the neighbourhood of 1°. If this were the 
case, one would either observe effects similar to those mentioned in the 
previous paragraph, or, if the relaxation time changed sharply with 
temperature an actual heat evolution when a relaxation time of a few 
minutes is reached. The absence of such effects seems to rule out this 
possibility as well. 

When de Klerk, Steenland and Gorter (1948, 1949) expressed doubts 
about the reliability of y-ray heating, they did so in order to account for a 
discrepancy between their specific heat values for chromium potassium 
alum below 0-1°K, and some of our preliminary results quoted by Hebb 
and Purcell (1937). Now, our early experiments were carried out on 
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cylindrical samples for which no 7* scale can be defined (Kurti and Simon 
1938 b) and the discrepancy between these two sets of results was therefore 
not significant. We nevertheless thought it worth while to carry out 
some new experiments on chromium potassium alum (Daniels and 
Kurti 1953) and although these showed none of the time effects or anom- 
alies mentioned above the results differed considerably from those of 
de Klerk, Steenland and Gorter (1949). This discrepancy is still un- 
accounted for and perhaps the most effective way of clearing it up would be 
to carry out successive experiments on the same specimens in the two 
laboratories. The use of chromium methylammine alum for such 
comparative experiments might even be more favourable as this salt has 
none of the undesirable properties of the potassium alum (Bleaney 1950). 

Finally, Professor Platzman mentions the chemical changes that are 
produced in the specimen by irradiation. This effect might indeed be 
serious when high y-ray intensities are used. In the calorimetric experi- 
ments we are concerned with, the rate of heat input is seldom higher than 
10 ergs/g/sec, and one ‘can estimate that even continuous heating for 
five hours would produce electronic excitation at the worst in only about 
1 in 10° of the polyatomic groups contained in the usual paramagnetic 
crystals. 

In conclusion it may be said that Professor Platzman’s detailed analysis 
is of great help in the assessing of the various possible errors that could 
be associated with y-ray heating. And, although present results show no 
effects that can be attributed to delayed heat evolution, experiments 
specifically designed for the purpose may well provide quantitative data 
on the mechanism of y-ray absorption in dielectric crystals. 
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SUMMARY 


In discussing the propagation of very low frequency radio waves to 
great distances, it is convenient to treat the space between the earth and 
the ionosphere as a wave-guide. This treatment is applied to some 
experimental measurements made by Dr. Weekes on signals of frequency 
16 ke/s, over a range of distances from 340 km to 3640 km from the 
sender. It is shown that the measurements can be explained in terms 
of the four least attenuated wave-guide modes. If the earth’s magnetic 
field is neglected, and the ionosphere is assumed to be a sharply bounded 
homogeneous, ionized medium, then it is shown that the height of the 
boundary is 69°1--0-5 km. 


$1. INTRODUCTION 


In the theoretical discussion of the propagation of very low frequency 
radio waves to great distances, it is often convenient to treat the space 
between the earth and the ionosphere as a wave-guide, and to consider 
the propagation in terms of a number of wave-guide modes. This was 
essentially the method followed in many well known theoretical papers, 
for example, Watson (1919), Eckersley (1932), Bremmer (1949). 

In the present paper the wave-guide treatment is applied to some 
recent experimental measurements made by Dr. Weekes, in which a 
receiver in an aeroplane was used to measure the strength of the signal 
from the sender GBR at Rugby (frequency 16 ke/s), over a range of 
distances from 340 km to 3640 km from the sender. It is shown that 
this treatment provides a simple interpretation of the results, and 
enables an estimate to be made of the height of reflection, and of some of 
the properties of the ionosphere. 


§2. Previous THEORETICAL WorK 


In two previous papers (Budden 1951 a, 1952), hereinafter referred to 
as I and I, the application of the wave-guide method to the propagation 
of a radio atmospheric was discussed, and it was shown that some of the 
experimental results (Gardner 1950, Bowe 1951) were conveniently 
interpreted by this method. In both I and II the ionosphere was 
assumed to be a sharply bounded homogeneous medium. In Paper I the 
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- relation between the wave-guide treatment and the alternative treatment 
in terms of successively reflected waves was discussed. It was shown 
that at short distances from the transmitter, the received signal is most 
conveniently considered as composed of a ground wave and a few 
reflected waves, but at great distances this picture becomes complicated 
because it is then necessary to include a large number of rays successively 
reflected at the ground and the ionosphere. If, however, the signal at 
a great distance is analysed into wave-guide modes, usually only a small 
number of modes of low order need be included because the higher order 
modes are heavily attenuated. In Paper II the full mathematical theory 
was given, and the effect of the earth’s magnetic field on the wave-guide 
modes was discussed. It was found to have little effect on the general 
form of the results, and in the present paper it is neglected. 

In Papers I and II the electro-magnetic field in a wave-guide mode 
was thought of as composed of two component plane waves, with their 
wave normals inclined at a (complex) angle 6 to the vertical. One 
component wave is converted into the other by reflection at the ground 
or at the ionosphere. After a part of a wave front has undergone two 
reflections, it is parallel to its original direction, and the phases and 
amplitudes of the original and of the twice reflected wave must be the 
same. This implies that only certain discrete values of 6 are possible. 

The changes of phase and amplitude of a wave can be represented as 
changes of a single complex quantity which is the complex amplitude of 
the wave. This changes both when the wave traverses a given path 
and when it is reflected at a surface. The earth’s surface is assumed 
to be a perfect reflector so that the only changes of complex amplitude 
arising from reflection occur at the ionosphere, which is assumed to 
have a (complex) reflection coefficient R. The wave front of a twice 
reflected component wave coincides with its original position after 
travelling a distance 2h cos 6, where h is the width of the wave-guide 
(height of the ionosphere). This introduces a difference of (complex) 
phase 2hp (cos @)/c between the original and the twice reflected wave, 
where p is the angular wave frequency and c is the velocity of light. The 
corresponding change of the complex amplitude is exp{—2ihp (cos 6)/c} 
and the product of this with R gives the total change of complex 
amplitude arising both from reflection and from the path traversed. 
This product must be unity for a self consistent mode and hence 


R exp {—2thp (cos )/c} =exp (27in), iaee a L 


where n is an integer, and is the phase difference measured in cycles 
between the original and the twice reflected wave. The reflection 
coefficient R is a function of the angle 6. The integer n may be positive, 
negative or zero, and determines the order of the mode, Equation (1) 
must be modified if allowance is made for the earth’s curvature. This 
modification is given in detail in Paper II and is used in §5 of the 


present paper. 
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The value of the reflection coefficient R depends on the model adopted 
for the ionosphere. The model to be used in this paper is the same as 
that in Paper I, in which the ionosphere was assumed to be a homo- 
geneous, sharply bounded, ionized medium containing N electrons per 
unit volume, making, on the average vy collisions per unit time. The 
effect of the earth’s magnetic field is neglected. At the level where very 
low frequency waves are reflected, v is of the order 107 sec}, and it is 
therefore assumed that v>p. Then if p is the (complex) refractive index 
of the medium 212 %9,/9 0 J, 2 
where Py=4rNe*/egmv, 

e=charge on the electron, 
m=mass of the electron, 
€y=electric permittivity of free space 
and the reflection coefficient R is given by 
R= {2 cos 0— +/(w?—sin? 6)}/{u? cos 0+ +/(u?—sim? #)}. . (3) 
A combination of (1) and (3) gives the equation 
{u2 cos 6—/(u?—sin? 6) }/{u? cos 0+ +/(u2—sin? 6) } 
—=exp{2thp (cos @)/e+-27in}, . . = « » «) (4) 
from which @ may be determined. 

The amplitude of the signal in a single wave-guide mode is given by 
an expression containing a Bessel function, but provided the distance 
from the transmitter to the receiver is greater than one or two wave- 
lengths it is permissible to use the asymptotic form for the Bessel 
function, and the signal at the ground is then proportional to 

p-1/*(sin 0)? exp{—tpp (sin O)/e}) >. 2 eto) 
(see Paper I, Appendix A), where p is distance measured horizontally 
from the transmitter. 


§3. ComposirE Wave-Guipr Mopss 

In the specific examples discussed in Paper I, it was found that the 
modes of orders 0 and —1 had very similar characteristics so that the 
fields in these two modes were practically indistinguishable even at great 
distances from the source. In the same way the characteristics of the 
modes of orders +1 and —2 were similar. This result applies also to 
the calculations of the present paper, and it is therefore convenient to 
treat the two modes of orders 0 and —1 together as a single composite 
mode, which will be called ‘mode A’. The field in this composite mode 
is given by (5) where sin @ has an effective value sin 0 a It is shown 
in $6 that sin 6, is given with sufficient accuracy by the average of the 
values of sin @ for the modes of orders 0 and —1. Similarly, the modes 
of orders +1 and —2 are treated together as a composite ‘mode B’. 
This has an effective value of sin 6, denoted by sin 0, which is a suitably 
weighted average of the values of sin 6 for the modes of orders —2 and +1. 
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§4, ANALYSIS OF THE EXPERIMENTAL RESULTS 


The study of the reflection of radio waves from the ionosphere by 
measurements of the signal at varying distances from the sender was 
first applied to very low frequency waves by Hollingworth (1926), and 
was used extensively by Best, Ratcliffe and Wilkes (1936), Budden, 
Ratcliffe and Wilkes (1939), and Weekes (1950). These observers made 
nearly all their measurements within a few hundred kilometres of the 
sender, and were concerned with the field resulting from the interaction 
of the ground wave and one predominant reflected wave. The e.m-f. 
of the ground wave is inversely proportional to the distance p from the 
sender, and they therefore found it convenient to display their results by 
plotting the product ‘ signal strength x p’ against p. By this method the 
field of an unattenuated ground wave alone would be plotted as a 
horizontal line, and the departures from this line could be ascribed to 
the interaction of a reflected wave with the ground wave. 

At great distances from the sender the field must be that of the least 
attenuated wave-guide modes, as explained in Paper I. The expression (5) 
shows that the amplitude in one composite mode is inversely propor- 
tional to y/p. This suggests that the product ‘signal strength x y/p ’ 
be plotted against p. For one composite mode the result should be a 
smooth curve showing an exponential attenuation arising from the 
imaginary part of sin @ in the exponential term of (5). Any departures 
from this smooth curve can be ascribed to the interaction of more than 
one composite mode. 

In fig. 1 this method of plotting is applied to some results obtained 
by Dr. Weekes from measurements made in an aeroplane which flew 
from England to Cairo and back in June and July 1950. The field 
strength from the sender GBR (16 ke/s) was measured continuously 
during the flight. Inspection of the curve suggests that there are two 
predominant wave-guide modes which interact to give an interference 
pattern. This suggestion is strongly supported by the fact that the 
spacing of the maxima and minima of the pattern is roughly constant, 
as it’ should be if the signal was produced by the interaction of two 
waves propagated radially outwards with constant but different phase 
velocities. The results displayed in fig. 1 were analysed on the 
assumption that the two composite modes A and B mentioned in the 
last section combine to give the observed signal. 

With this assumption it is possible to obtain three numerical 
quantities from the results shown in fig. 1. These are: 

(i) The difference of the reciprocals of the phase velocities V, and Vy, 
for the modes A and B. This gives a value for # (sin @,)—A(sin Gy), 
where the symbol # denotes ‘ real part of’. 

(ii) The attenuation constant of mode A. This gives -% (sin @,) 
where .% denotes ‘imaginary part of’. 

(iii) The attenuation constant of mode B. ‘This gives SF (sin 6p). 
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The relation (i) was found as follows. Successive maxima and minima 
in fig. 1 correspond to increments of 7 radians in the phase rahe 
between the two modes A and B. The positions of the maxima an 
minima were therefore plotted as abscissae with ordinates increasing in 
steps of 7 radians. The resulting graph is a fairly good straight line, and 
its slope in radians per unit distance is easily shown to pg equal to 
p(1/V,—1/Vy). Now from (5) 1/V,= (sin Oa)/e; 1/Vzp=2 (sin Pale, 
so that the slope of the line is equal to {# (sin@,)—F (sin Oy) }/c. It 
was found that & (sin 0,)—Z (sin 0g) =0-0181+0-0003 (standard error 
estimated roughly by inspection of the graph). 


Fig. 1 


Signal amplitude X 4/(distance)—> 


~ 1000 8.8 S3N00 3000 


Distance in km > 


The product ‘signal strength square root of distance’ is plotted against 
‘ distance ’ for the signal from the sender GBR (16 ke/s) during the day- 
time in June and July, 1950. 


The attenuation constant of mode A ((ii) above), was found by 
smoothing the curve of fig. 1, to eliminate the effect of mode B. Let 
E,, LE, be the ordinates of two adjacent maxima (or minima) in fig. 1, 
and let H, be the ordinate of the intervening minimum (or maximum). 
Let H,=}(H,+#;+2E,). Then £, was taken as the smoothed ordinate 
for the distance corresponding to E,. This is a somewhat crude method 
of smoothing, but the accuracy of the experimental data is not high 
enough to justify a more elaborate method. Log, H, was plotted against 
distance, and a straight line drawn through the points. The slope, xq, 
was taken as the attenuation constant of mode A. From (5) it is clear 
that y,=¥ (sin @,)p/c. It was found that / (sin 64 )=0-00048+0-00003 
(rough estimate of standard error). Using this value of .% (sin 04), the 


Y 
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quantity “signal amplitude x /p’ in mode A alone was found, and is 
shown as a dotted curve in fig. 1. 

The departures of the maxima and minima from the dotted curve in 
fig. 1 give an indication of the quantity #,—‘ signal amplitude x \/ ps 
for mode B alone. Log, Hy, was plotted against distance, but the 
resulting points were so widely scattered that only a very rough estimate 
of the attenuation constant y, of mode B could be made. yg can algo 
be estimated, when x, is known, by comparing the signal amplitudes 
in modes A and B at the same distance from the transmitter, because 
the theory shows that the two modes should have approximately equal 
amplitudes when p=0. This method also gave only very approximate 
values of yz, but the two estimates were consistent. It was found that 
F (sin 6,3) probably lies within the range 0-002 to 0-003. 


§5. SUMMARY OF THE THEORY OF Wave-GuiIpE MopEs 


It is now desirable to find a model for the ionosphere, in which 
- composite modes A and B would be propagated with the characteristics 
found in the last section. It was explained in §2 and in Papers I and II 
that some progress with the theory is possible, if the ionosphere is 
assumed to be a homogeneous sharply bounded medium, and it will be 
shown that this assumption leads to a consistent explanation of the 
results. 

First, eqns. (1) to (4) were solved to give the characteristics of the 
simple modes of orders —2, —1, 0 and +1. It is convenient to put: 


Dia ip je—pok=e *; cos0=Co 2... (6) 
Then eqns. (2) and (3) give 
er {(y4—1)C— »/ (uC? iu) }/ {(u—1) 0+ »/ (we? —tu)}. . (7) 
If the earth is assumed to be flat, eqn. (1) gives 
oes ar te ewe ov 8) 


If the curvature of the earth is allowed for, and K denotes its effective 
value, then it was shown in Paper II (loc. cit. eqn. (61)) that provided hk <1: 

a—n7=BuC—tKhpu/C. 2S a aoa ()) 
(7) and (9) are two equations which determine C’ for given values of 
B, w and n. They were solved with the aid of the EDSAC calculating 
machine in the Mathematical Laboratory of the University of Cambridge. 
Details of the method are given in the next section. 


§6. MeTHOD OF CALCULATION 
Equation (9) is a quadratic for C and the required solution is 
C=(a—nr)/2Bu+t {(a—n7)?+2KhB wu? }47/2Bu.  . . (10) 
Values were first chosen for 8, wu and n. A trial value of « was then 
taken, and (10) used to find C. This was put into (7) to give a new 
trial value of «. The process was then repeated, and was found to 
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converge. The EDSAC programme was so arranged that the machine 
compared two successive values of a, and when the difference had 
become less than a specified quantity, the calculation was stopped, and 
the values of u, 8, n, «, C and sin@ were recorded. The calculation 
was then repeated with a new value of wu. From the results it was possible 
to tabulate F (sin @) and .% (sin @) for a series of values of u, 8, and n. 

Let 0), 6,, be the calculated values of @ for the modes of orders 0 
and —1 respectively. Then the signal in the composite mode A is 
proportional to 

2S ,—p~1/2{(sin 6,)!/? exp (—ipp sin 4,/C) 
+-(sin 0_,)1/? exp (—tpp sin@_,/C)}. . . . (11) 
Let sin 6)-sin 6_,=2 sin 0, ; sin@,—sin@_5=24,, ~. = (12) 
where 4, is small. Then approximately from (11) 
S,—=p1/? (sin 0,)1/? exp (—ipp sin 0,4/C)[cos (pp4,/C) 
| —iiA,(sin 6,)-1 sin (pp4,/C)]. 2, Give aa oe ies 
From the numerical results for the largest value of p used, it is found that 
| cos (pp 4,/C)| 0-97; arg {cos (pp44/C)} = 50’ 

and | $4, (sin 0,)~! sin (pp 4,/C) | 0-00024. 
The second term in the square bracket of (13) was therefore neglected, 
and the first was set equal to 1. 

For the composite mode B a more complicated process was necessary, 
because the difference of the attenuations of the modes of order +1 
and —2 was appreciable, and at 3640 km from the source the amplitude 


of the mode of order —2 was only about one-half that of the mode of 
order +1. The signal in the composite mode B is proportional to 


2S},=p/? {(sin 6, ,)!/2 exp (—ip sin 6, ,p/C) 
+ (sin @_ ,)!2 exp (—1ip sin @_,p/C)}. . . . (14) 
The factors (sin@,,)!/ and (sin @_,)!/2 were both sufficiently close to 
unity to treat them as equal. It was only in the exponentials that small 
differences between sin @,, and sin@., were important. Let 
exp (—zp sin 6, ,p/C)+-exp (—ip sin 6_,p/C)=exp (—ip sin 0, p/C). . (15) 


Sin Or could be calculated from (15) for each value of p but it was found 
sufficiently accurate to take 


ee 
2sin@,;=asin6,,+bsin@., .. . . . (16) 


where a+b=2, and a and b are constant weighting factors, which were 
calculated from (15) and (16) for p=2000 km. This is the value of 
about half way between the extreme values shown in tgs The 
resulting error in phase for the largest value of p used was then about 20° 
and an error of only 67 km in assessing the position of the last ain 
in fig. 1 would give a similar error in phase. 
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For a given value of f, # (sin 0,)—Z (sin 0,) was plotted against w, 
and the value uw) of w was found which made the ordinate equal to the 
experimental value 0-0181. For this value of wu, % (sin @,) was then found. 
The process was repeated for other values of 8, and the results are shown 
im fig. 2. The experimental value of .% (sin @,) was 0-00048, and it is 
clear from fig. 2 that this is given by B=10. In this way it was possible 
to find the values of 8 and u which fitted the experimental results best. 
The appropriate value of uw) was found to be 2-35-40-02, and since 
p=2r X16 rad. sec”! this gives p,=4:28x104 rad. sec-?. With these 
values, the second relation in (11) shows that h=69-1-L0-5 km. For the 
above values of 8 and w, the value of % (sin 6,) was found to be 0-0027, 
which is within the range indicated by the experiments. 


Fig. 2 
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y (curved earth) 
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0-0004. 


JF (sin 6,)—> 


Theoretical Value 
(flat earth) 


0:0002 
4 


Value of B > 


Determination of 8 from experimental results. The experiments show that 
F (sin 04)=0-00048. 


§7. Discussion or RESULTS 


It is natural to ask whether the results of the previous section would 
be appreciably altered if the curvature of the earth were not allowed for. 
Some calculations were made in which the earth was assumed to be flat, 
and the results are shown in the lower curve in fig. 2. This is seen to be 
consistently below the experimental value for / (sin @,), and no reasonable 
values of 8 and w could be found which would fit the experimental 
results. It is therefore essential to include the effects of the earth’s 
curvature. 

The value 69-1+0-5km found for h, the height of the ionosphere, 
depends only on the experimental value of # (sin 0,)—& (sin @,) and 
not on ¥(sin@,). This is because the product Pu) proves to be 
independent of f for a range of values of 6 around 10. #(sin 6,)—Z (sin 0) 
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is the most accurately determined of the three numerical quantities 
deduced from the measurements, and hence h is given within narrow limits. 
It is interesting to compare this value with that given by Bracewell et al. 
(1951) who deduced a height of 74km from measurements made on 
summer days at steep angles of incidence. Since the component waves 
in the wave-guide modes impinge at nearly grazing incidence on the 
ionosphere, the two values obtained for the height of reflection would 
not be expected to agree exactly. 

In § 2 an expression for p, was given which shows that it is proportional 
to the ratio of the electron density N to the collision frequency v. The 
best estimate of v at a height of 70 km seems to lie between 1 x 107 and 
3x 107 sec“. Then, using the value of p, given in the last section, it is 
found that N lies between 135 and 400 cm-?. 

It is of interest to find how far the fields of the wave-guide modes 
penetrate into the homogeneous medium which is assumed for the 
ionosphere. The upgoing component wave of a mode is incident on 
the ionosphere at a (complex) angle 0, and gives rise to a transmitted 
wave in the medium with its wave normal at a (complex) angle ¢ to the 
vertical. The vertical variation of the field in the medium is thus given by 
exp {ip (t—cos ¢ . wz/c)} where ¢ is time, and z is height. The amplitude 
in the medium is therefore proportional to exp {—p|.4(p cos 4) |z/c}. 
It falls to half value in a vertical distance A (In 2)/27.4(y cos d), where A is 
the wavelength in free space. The value of this distance was about 
4-7 km for all the four modes of orders —2, —1, 0 and +1. Thus it is 
probable that only a thin region of the ionosphere plays an effective part 
in the propagation of low frequency waves to great distances. 

Figure 1 shows clearly that the oscillations of the curve at small 
distances from the transmitter are less smooth and regular than at 
greater distances. This could be explained by the presence of modes 
of higher order. These would be more heavily attenuated than the 
composite modes A and B, and so would cease to have appreciable 
amplitude at greater distances. 


§ 8. CONCLUSION 

The propagation of very low frequency radio waves (16 ke/s) to 
distances between about 400km and 4000km can. be adequately 
explained in terms of wave-guide modes, and to account for most of 
the experimental results only two composite modes (four simple modes) 
need be considered. If it is assumed that the wave-guide is bounded 
by a perfectly conducting ground, and an ionosphere which is a sharply 
bounded homogeneous ionized medium, then good agreement with the 
experimental results can be obtained, and estimates can be made of the 
height of the ionosphere, and of the electron density in it. It is found 
that the height is 69-1-+.0-5 km and if the collision frequency is taken to 


be between 1107 and 3x 107 sec-! at that height then the electron 
density is between 135 and 400 em-3, 
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ABSTRACT 


Interferometric and optical studies on a twinned diamond reveal the 
presence of true crystallographic slip on two opposite faces of the twin. 
Slip occurs parallel to a (111) plane, in both cases, making thus an 
angle of 704° with the surface. An internal optical opacity is found 
associated with the slip. The two faces on either side of one slip line are 
not parallel but inclined at the very small angle of 10~° radian. 


§ 1. INTRODUCTION 


Ir is doubtful whether earlier observers have been able fully to 
substantiate the existence of slip in diamond. Williams (1932) has 
suggested that certain features appearing on some microphotographs of 
diamond might be attributable to slip, but the evidence is somewhat 
tenuous. Since diamond cleaves only effectively on the (111) octahedron 
plane, it might be anticipated that if slip were to be found, then it should 
occur parallel to this natural cleavage plane. 

During an extensive interferometric study of growth features on the 
surfaces of many diamonds, we have found on one crystal, a contact 
twin (i.e. ‘ macle ’) loaned to us by Dr. Grenville-Wells, decisive evidence 
of crystallographic slip which we have studied quantitatively by 
interferometric and light profile methods (Tolansky 1952a,b). Slip 
regions have been found on both the major (111) faces of the twin and 
it has been established also that the slip plane is indeed parallel to (111). 

To study this crystal use has been made of a number of optical 
techniques. The surface had on it some high projecting regions and 
since it has long been established that multiple beam interferometry 
requires close approach to the two surfaces involved, a special micro-flat 
technique was developed, in which the flat used was the small tip of a 
truncated cone, which could be brought into close approach with the 
desired region. This technique is being described elsewhere. Again, 
high magnification multiple beam interferometry was desirable in some 
regions and we reproduce later an interferogram taken at x 2000, 
using for this the thin film technique we have already described 
(Tolansky and Omar 1952) elsewhere. 

a Senne ne eee eas neem SE 
*Communicated by the Authors, 
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§2. OBSERVATIONS 
(a) Slip on the First Face 


On the one side of the crystal the slip appears as a straight line, cutting 
right through all the growth features, some 15 mm long, running right 
across the surface. Part of this is shown in fig. 1, Plate 29 (x 30). The 
series of parallel lines KL are growth sheet fronts, and at 60° to these 
lies XY the slip line. The topography of this region is rendered clear 
by fig. 2, Plate 29 (x 100) which is a multiple beam interferogram. 
Here MN is a growth sheet step, parallel to the group of steps KL and 
once more XY is the slip direction. From the fringes one can readily 
evaluate the topography, a simplified diagram being shown schematically 
in fig. 7. The Fizeau fringes of fig. 2, Plate 29 are ambiguous as to step 
directions and these have been clarified by using fringes of equal 
chromatic order (Tolansky 1945). With these fringes it is confirmed 
that in fact the level differences A-~B and C—D are identical in the 
neighbourhood of the cross over. 


Fig. 7 


EG 
Sy 


An examination of the slip step over the entire length of the crystal 
reveals the existence of a slight slope. The step increases regularly 
from the left to the right. Thus at the extreme left-hand side the slip 
is only 480 a. At a point 2mm to the right it has increased to 680 4, 
and at 8 mm to the right to 1280 a. This corresponds to a slope between 
the slipped regions with the very small angle of 10~° radian, i.e. merely 
2 seconds of arc. Although small, this value is precise and has been 
confirmed by measurements on intermediate positions between ends. 
If this slip were to be taken as due to a dislocation the angle would 
imply a displacement of one atomic lattice in 10° lattices, to account 
for the slope. However, it should be pointed out that the slip does not 
run to zero at one end of the crystal, the minimum step of 480 a clearly 
corresponding to a slip of some hundreds of lattices. 

Examination of the crystal in transmission with polarized light reveals 
a general distribution of strain other than in the slip region. Near the 
slip region itself and undoubtedly closely associated with it is a peculiar 
opacity. This is shown in fig. 4, Plate 29 (x 90). Examination with an 
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intense light source shows that the opaque region is sharply defined by 
the slip edge and then falls off somewhat, the final edge (lower) tailing 
away in a fuzzy manner. Inside the opaque region are some faintly 
transparent interference fringe regions. The difference in sharpness at 
the edges of the opaque region is clearly shown in the picture. 

It will be shown later that the slip makes a large angle with the 
surface. This being so, high magnification experiments were undertaken 
to see whether the slip step discontinuity was sharp. To do this a 
silvered thin film technique (Tolansky and Omar 1952) was used. 
The crystal was silvered for multiple beam interferometry. A very thin 
film of Canada balsam was drawn over the slip boundary by means of a 
dilute solution of xylol. This dried rapidly and was then silvered. 
Multiple beam fringes are then formed in the thin film, the optical step 
being, of course, magnified by the ratio of the refractive index of the 
film to that of air. We have already demonstrated that high lateral 
magnification multiple beam fringes are obtainable by this procedure, 
fig. 3, Plate 29 shows the fringes across the slip step at x 2000. When 
the high lateral magnification is taken into account, it will be recognized 
that the fringe definition is exceptionally good. 

Apart from the fact that such a picture permits accurate determination 
of the step, with an error of but a few angstrom units, the picture 
establishes that the step is discontinuous, sharply so. Now evidence 
will be produced later that the slip steps found on both faces have an 
undercut angle of 703° and there is every reason to believe, that the same 
will always be true. The sharp discontinuity observed is then only 
what can be expected if indeed the slip step undercuts, and is thus far 
a satisfactory observation. 

It is to be noticed that the slip line is exactly at 60° to the growth 
direction and is thus parallel to an edge of the (111) face. 


(b) Slip on the Second Face 


Because of the twinned nature of the crystal the slip cannot be 
expected to penetrate from one side to the other and indeed this is the 
case, for the line XY does not appear on the opposite face of the flat 
stone. However an independent slip line does appear on the second 
face, this being inclined at 60° to the direction of XY. This second slip 
line also traverses right across the crystal and in doing so very clearly 
reveals a feature of some considerable interest. For the second face 
has on it a number of high growth plateaus, as well as deep triangular 
depressions (familiar features on (111) faces). The slip in question runs 
across some of these and in doing so becomes quite clearly displaced. 
It is from these displacements that the angle of undercutting can be 
determined, as shown in fig. 8. Suppose we have a growth step P. If 
this be cut by a slip plane normal to the surfaces then the cut lines in 
the upper and lower growth regions, when viewed from above, are 
colinear, If the cutting slip plane has an angle to the surface, then the 
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two lines will appear displaced, as in fact shown schematically in fig. 8. 
If the growth sheet step height be known, and is reasonably large, the 
measured displacement in the slip line permits the angle of undercut 
to be calculated. Furthermore the direction of displacement determines 
the direction of undercut. An example of a measurement made this 
way is shown by fig. 5, Plate 29. Here the slip traverses a growth trigon. 
The discontinuity in the slip line as it passes over the growth step PQ, 
for example, is readily seen. This step height PQ measured by the light 
profile is large, i.e. 90. The lateral displacement of the slip line 
measured in fig. 5, Plate 29 (x 30) is 32. From a simple calculation 
it follows that the angle of undercut is 704°. This is precisely the angle 
made by adjoining (111) faces on diamond. We have therefore 
established that the slip line is not only parallel to a (111) edge, but 
aa the whole slip plane is parallel to a (111) plane, i.e. to a cleavage 
plane. 


Fig. 8 


This interesting and exact geometrical deduction has only been made 
possible because of the occurrence of one very large discontinuity (90 1) 
over which the slip traverses on the second face. Only in one position 
does the other slip line on the first face also pass over a relatively large 
discontinuity, and it is striking that from this face also it has been 
established to have an angle of 703°, exactly as in the case of the second 
face. The light profile picture over this part is shown in fig. 6, Plate 29 
(x 90) in which the growth discontinuity is 31-5 and the slip line 
shift lly. In fig. 6, the slip line is that which runs in the north-east to 
south-west direction crossing the group of parallel growth steps, the 
discontinuity appearing markedly in the last step. The optical light 
profile, which also shows the discontinuity at this step, is the vertical 
line having on it a short horizontal cross, some 1 cm from the top of the 
picture. The slip on the second face has also associated with it an 
internal opacity. Furthermore the second-face slip dose not penetrate 
the two twins but, as in the first case, ends within the body of the crystal. 


SER. 7, VOL. 44, NO. 352.—MAY 1953 2M 


518 On the Observations on Slip Found in a Diamond 


(c) A Possible Explanation of the Opacity ; 


Now that it is established that the slip plane is making an angle 
of 704° with the surface it is quite possible to explain the opacity on the 
following lines. If the slip has developed partly into a crack, then an 
air film in the crack would totally reflect the incident light at the 
interface of diamond-air. The incident light would make an angle of 703° 
with the perpendicular to the slip plane which is greater than the critical 
angle of diamond (243°). As either part of the macle is 1-5 mm thick 
and the opacity at its maximum as estimated in different parts is } mm, 
the crack, on this assumption, could only have penetrated to no more 
than two-thirds of the thickness of either part. If this interpretation 
of the opacity is correct it is reasonable to assume that the crack might 
have a finite width at its mouth where it emerges on the crystal face. 
Although the formation of such a small wedge crack might lead to the 
appearance of a small inter-facial angle between the two parts on either 
side of the crack, no such very small inter-facial angle has been detected. 
It is clear that the observed step cannot be accounted for by postulating 
such a mechanism for this suggested alternative would lead not to fringe 
steps but to fringes showing an angle kink. The fringes would appear 
continuous but change abrupt direction at the crack and would not, as 
are observed, be violently out of step. . 

It is emphasized that the observed angle of 10-° radian referred to 
above applies to the fact that the slip is emerging out of the surface 
with a slightly increased step height as one traverses along the length 
of the ledge. 

Thus it may be concluded that 

(a) the occurrence of crystallographic slip on diamond has been 
demonstrated on two faces of a twin, 

(b) the slips takes place on (111) planes, which are the planes of easy 
cleavage, 

(c) the two faces on either side of one slip are inclined at the slight 
angle of 10~-° radian, 

(2) the slips have associated with them regions opaque to transmitted 
light, 

(¢) the slip appears on both sides of a twinned macle, but in neither 
case does the slip penetrate to the opposite side. 
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ABSTRACT 


The vibrations of a bounded lattice are identified with those of an 
elastic medium (isotropic, arbitrary elastic constants) having stress free 
surfaces. The densities of the representative points in wave number 
space for the four possible types of vibration are estimated. These 
estimates are carried one step further than the usual (Debye) method by 
including terms proportional to the surface. Summing over the densities, 
the total number of excited modes and the intrinsic energy of the lattice 
vibrations are obtained in the usual manner. From these we calculate 
surface terms for the Debye temperature and specific heat and for the 
free energy of the lattice vibrations. 


§1. INTRODUCTION 


THE specific heat of a solid has been calculated (Debye 1912) by replacing 
the thermal vibrations of the lattice with elastic waves in a continuous 
and isotropic elastic medium of the same volume. For a large volume, 
arbitrary surface conditions are selected and the total number of modes 
M(w,,) having a frequency (v,,=a,,/27) calculated. w,, is then deter- 
mined by equating M(w,,) to the total number of degrees of freedom of 
the lattice points (ions, atoms or molecules), ic. 3NV, where N is the 
number of lattice points per unit volume and V is the lattice volume. 
If, however, the ratio of surface area to volume becomes large, natural 
surface conditions should be used thus affecting the value of w,, which 
becomes surface dependent. The natural condition appropriate for a 
continuous medium will be a stress free surface. 

The intrinsic energy, which corresponds to the sum of the average 
energy of the vibrations over the M(W,,,) modes, thus contains a surface 
term from which a surface term in the specific heat and a temperature 
dependent surface free energy can be calculated. 

The calculation has been carried out by Brager and Schuchowitzky 
(1946, referred to as B) for the case of an incompressible medium. This 
gives a great simplification since only transverse waves can occur. They 
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considered the two-dimensional vibrations of a parallel plate (normal to 
the plane of the plate) of infinite extent but finite thickness and then 
extended the solution to three-dimensional waves. Their treatment is 
however in need of correction since they pay no special attention to. 
waves with a zero wave number component in the direction of the thick- 
ness of the plate and carry out the extension to the three-dimensional 
case incorrectly. 

Besides investigating the effect of these corrections, it seems desirable 
to extend the calculation to the case of arbitrary compressibility, thus 
including longitudinal waves. This not only gives a better insight into 
the present problem, but consideration of the longitudinal waves is 
required for the calculation of the electron-lattice interaction in a finite 
plate. This last problem arises in the discussion of the surface energy of 
a super-conducting specimen to which we hope to return in a future paper. 

Our calculation is based on the representation of the possible modes of 
vibration, for the various types that are found to occur, by points in wave 
number planes. The densities of these points are then estimated and the 
total number of modes obtained by summing over the densities. Having 
determined the Debye frequency we can calculate the intrinsic energy, 
again using the density of modes. 

In §2 we consider two-dimensional waves first as they are simpler to 
interpret and occur as special cases of the three-dimensional waves 
subsequently treated in §3. The evaluation of the surface terms for 
the specific heat and free energy of the lattice depends only on the final 
result of § 3 which expresses the total number of excited modes of vibration 
in terms of the maximum frequency and elastic parameters. 


§2. THE Two-DIMENSIONAL CASE 
We consider a parallel plate of infinite extent but finite thickness which 
fills the space between z= +L, and z=—L, of a rectangular co-ordinate 
system. Solving the elastic stress-strain equations, we can express the 
displacement in the (x, z) plane as (cf. Lamb 1917) 


where the potential functions ¢ and % are given by 


p= — USF, cos q,,2 exp [i(q,x—at)] | ‘cal 
TASS ee ree s symmetrical waves ede: 
p= 2°P, sin g,,2 exp [1(¢,7%—at)] (2) 
or 
9A, a Bs s . ’ 
p= 2*Ff sin 7/2 exp [i(q,7— at) | : ' 
les antisymmetrical waves 3 
p= —2i*P, cos 4,2 exp [1(¢,7%—wt)] Fs a (3) 
with respect to the ay plane. 
P : : 
"he wave number components satisfy the relationships 


Ge t+Ie=G7 = w/e? = ww po/ (A+ 2h) 
Ge Ia =9f S w/e? = co pol. } 
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Here w is the angular frequency of the elastic waves whose velocity 
(c;—longitudinal, c;—transverse) is given in terms of the usual elastic 
constants (wu and A) and the density p, of the material. The form of the 
constant coefficients will be shown to be convenient in later work. 

In the x direction we take the usual periodic boundary conditions with 
periodicity 2L,(L,>L,) so that 


NT ‘ : 2 
Umea free! (USM 3 4, Ve, OE ee Ae ae ee)s os (5) 
x 
while the surfaces z= -+ L, must be stress free. The latter condition gives 


the spectrum of the waves; namely 
tan@, 4(1+ K)a(1+2K—a?K) | 6 
tan Q, (12K +02)? ee ey ee 
for the symmetric and 


tan@, 1 (7) 
tan Q;, USCA P CN | Et eh ee 


for the antisymmetric waves. Here we have introduced the quantities 


QO:=Iel: > Q2=Ubh, : %=Q (Vi > K=p/a, Seer ee (8) 
where K can be expressed in terms of Poisson’s ratio 
K=p/A=(1/20)—1. Sa}. Bo, tie Aa Semen Ge, 
Further from eqns. (4) 
OR 2 OU) ESV OOS AGE Gs 0 0 < . (10) 
WOO KU 2 kK og I+ KD. Ce CS) 
where 
OF eM ety ie eee EZ) 


‘For any pair of solutions Q,, and Q,, satisfying (6) or (7) the ratio of 
the coefficients will be given by 
(AP) sin 2Q), 
OUR RE os tf eae 
(GP, sin 2), 
There are three possible types of wave motion : 
(1) ‘ Pure Rayleigh Waves ’ 


(13) 


de de -- Viz and gq, are imaginary. 
(2) ‘ Mixed Waves ’ 
. GP >Ue> Ue. Wz is imaginary and q, is real. 
(3) ‘ Pure Bulk Waves ’ 
Gd. + Ty and q, are real. 


‘From the displacement function it will be seen that waves of type (1) 
are confined to the surface of the plate, those of type (2) consist of the 
superposition of a longitudinal surface wave and a transverse bulk wave, 
whose amplitude is constant over the volume of the plate, while those of 
type (3) consist of the superposition of longitudinal and transverse bulk 
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waves. We now discuss the frequency spectrum for the three types and 
determine the number of modes of vibration having a frequency less than 
a fixed maximum frequency (the Debye frequency). 


2.1. Pure Rayleigh Waves 
Introducing the real quantities 
0 .=-—9" and Q' =—92", : e ° sd (14) 
so that c=Q,./Q,.=' ./Q’ . (taking the positive square root), the spectrum 
for symmetrical waves becomes 
tanh (ii Lon OC eee ae Lae 
tanh Q’,, tanh Q’;, 
First of all there exists a solution of (15) as Q’;,->00, namely the real, 
positive solution of the cubic 
Ix(%)+1=0, fo 
say %)(K). (There is a factor («—1) for all K which corresponds to zero 
frequency.) It can be easily shown that 
aq(0) <a9( A) <ao( 00) 
[ax9(0)=0-2956, a9( 00) =0-6181]. 
a 9(A) increasing monotonically from one extreme to the other. 
Turning again to eqn. (15) we consider the solution for « > 0. This is 
Wn _ 41+K) . 
tanh Q@’,, 142K ’ 
which is approximately Q’,, & 4(1+K)/(1+2K). Finally, for large 
Q', the asymptotic solution is 


a=0 iLeQ'=0, 7 See 


2 
—>—— exp (— 2« "|. nce Rees 
g (a) i ( o& lz ( ) 
(97’x(%) is always negative]. If we now exhibit the solutions of (15) as 
a curve in the Q’,., Q’;, plane the possible frequencies will be given by 
the intersections between it and the family of hyperbola defined by 
eqns. (10) and (5). From (11), the maximum frequency w,, corresponds 


to a maximum @";, while (10) gives the range of q, in terms of that of 
Q’;,.. Thus the total number of modes is 


Q' Ost Qi: Oy 


Hsien) 1dy,o2 = oe (Oe 
7 2rL, 1+-K (1-+2K)(1—a2) 
* 1 (es 1/2 4(1+-K) . 
Imtlz \ 14K 1+2K |’ (19) 


where ¢,,;=w,,/¢, is the maximum transverse wave number and 4A =4L a 
nN aetor 2 Cd—rSs ] it] rae 
The ae occurs since Q,, can have negative and positive values. 
Now M;. is proportional to A/L., i.e. the surface area (for two-dimen- 
sional motion A represents the volume and A/L, the surface area, if we 
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consider unit thickness in the y direction). Thus we can neglect the 
second term in the bracket which corresponds to the non-zero value of 
Q’,, given by (17). 


A 
M;S(o,)= QrL, YoGint ; Yo= 


A similar argument applied to the antisymmetric waves, where the only 
difference is the lower limits of the q,, range, gives 


M4 (w n)= M,§ (Wm)> 
, A 
so that Mw n= M4 (w m+ M8 (Om) = ae Ps Yodint: Go Big (21) 
aL, 
2.2. Mixed Waves 


Since q;, is imaginary we introduce «’=i« so that the spectrum for 
symmetrical waves becomes 


tan Q v4 ay 
and". = — 1 p(—ta') =ape(a’), . 2 . . (22) 
2] 
where 
no 4(L Ko (12K Ko’) ie 
P(x’ )= eer eoKk=o2 he. (23) 
For large Q’;, the solution of (22) is 

Q,,=Npm/2+8; Stan! da’) Jo<o< 5 | Maha 754) 


injg= 02, 4, 6) 8 >. 
The possible frequencies will be given by the intersections of the curves 
represented by eqn. (24) and the ellipses represented by eqn. (10) with 
Q' 2=—Q/2. The total number of modes will be the number of inter- 
sections, inside the quarter circle given by eqn. (11) with w=w,,. The 


radius of this circle is 

a [poll K) i 
Pra WV K (142K) Sc ol a eee, Se 
To calculate the density (p) of the representative points we consider 
the intersections A, B, A’, B’, for consecutive values of n,, and n, 
(cf. fig. 1). The density will be given by the reciprocal of the area 
ABB’A’. In calculating this area we can neglect the curvature of AB 

since it introduces higher order approximations. 


Areas b bA= AB ERP=—DC.GF. . 7... (26) 
Now 
=) (% :) 7 . 
DC=|-—=) 2 and FG= meen TS ti, 
Cn Ux Od, an L ( ) 


Substituting into eqn. (26) the density is given by (factor 2 to account for 
positive and negative ¢,,) 


eal ania (ons OF (28) 
Pe DO EG 7 \ce Oy \O0'tc) Oe bie Os 
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Evaluating the derivatives ((24) and (10)) 


4a penis (ee 
P= oA AKI K)A|” 1K ae 1+2K{ |° 


The total number of modes is thus 


2 r-Pm 
MySon)=| |" eqdadg+% . - - + (30) 


where we have introduced the polar co-ordinates 
ve el 2 MA A — / i gn oe / . 
q =r +Q lz » tan d=Q,./@ Ee Pera (hic! 
and Z, is the correction for modes with n,,=0 which are included in the 


integral only with weight $, i.e. 


a (Pm D Ldn ‘ 
y= 3). (P)ayg-0 42! p= ital Be ieee Soak ee 


(This term is omitted in B.) 


Fig. | 


D C 
Q+tZ 


Neighbouring representative points for consecutive values of mz and ny. 


Substituting into (30) for p 


4 A> Frq* L+-Kk\1/2 
M y°(o,,)= ane | oe tan- 1} ( w =) — Imi (5 —n1(K)) | . . (33) 


K L, 


& 


where 


- lf -K: \V2 (928! (tan 6) /2K-+1 
ihe / =—=_ ————— K a o >) € 
o(K) © (; z) a ———_— ( K +tan? 6) dp. -. . (34) 


cos 4 


4p aay it > ; » ] You i . . . A 
. The value of the integral 7,(0) can be deduced from a similar integral 
given in B. 


me Uriel Pity cake Paar main, 61-9 
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We have not evaluated the integral for any other value of K since it does 
not occur in the more important three-dimensional case. 
The spectrum for antisymmetric waves becomes 


tan Q,,/tanh Q’,=—Ifb(a’)  . . . .). (36) 
so that again for large Q’,, 
Q c= NyT/2+6; S=tan-1dh («’) [(0<d8<n/2] . . . (837) 


mod i= 15.30.74. 

Proceeding as for the symmetric case the density of representative 
points is the same so that the only difference is the absence of the zero wave 
number correction 


A A Gn ; 1 =e ie Qmt™ = 

: el Gi a3 [ = fan 4 (=) a is THK )| tee oS) 
and adding 

My(w,,)= (Wm) +M x8 (Wyn) 


Cas 1-2K ue Vint™ 1 7 ‘ 
fi tar (Se ) + 8 (714K) |. ee eS) 


The approximation tanh Q’;.~1, which has been used throughout, 
implies Q;,>1. This is satisfied except for a narrow region of area P,, 
near the Q’;, axis. The value of p given by (29) tends to zero with Q’;, 
whereas the actual value is of order A/L,’._ Thus the relative number of 
states not satisfying the condition is of order P,, A/L?M(w,,)2O(1/q,,L- 
Since the approximate solutions are themselves used to calculate a surface 
effect (i.e. 5 occurs only in the surface terms) the approximation is justified. 


2.3. Pure Bulk Waves 
The formal solution of (6) may be written in the two ways 
Ooo; 0, tana |g tan@,|- .;. « . (40) 
or pee 07 who == tans { (lig) tan Oy) me 4 (41) 
where the n are non-negative integers and 0<6, ;<7. This case differs 
essentially from the one previously treated due to the fact that the 6 
contain rapidly varying trigonometric functions. 

We first of all consider the representative points in the Q,,, Q;, plane. 
The solutions of (6) giving the relation between Q,, and Q,, are represented 
by the curved lines passing through the checkered pattern of fig. 2. For 
real values of Q,, (cf. (10)) 

a? > (14+2K)/K o. g(a) SO. 
Thus there are no solutions inside the shaded squares. The representative 
points are given by the intersections of these curves with the family of 
hyperbola (10). The total number of modes will be the number of repre- 
sentative points inside the figure formed by the straight lines «=0 and 
x=(1+2K)'/?/K1/, and the rectangular hyperbola defined by (11) with 
Wa= Win: 
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Proceeding as for type 2, the density of representative points can be 
expressed in the two forms 


cz ~29(55) (4) eae geen ees 
Pr \8Q a) on \OOte) 042 
L,, {on 0g : 
s_ dx ( Ne zx he eS 
ae Cae Cae 
depending on whether we project AB of fig. 1 onto the Q;, or Q ), AXis. 
Fig. 2 


B is one of the family of hyperbola defined by (10). The curved lines passing 
through the blank squares of the checkered pattern represent solutions of 
(6) for symmetric (full lines) and (7) for antisymmetric waves (dotted 
lines) respectively. The intersections of B with these curves are 
representative points for waves of type 3. A is the bounding rectangular 
hyperbola given by (11) corresponding to the maximum frequency. 


Considering p,* first and substituting from (40) 


A K \12 /1+2K 1+-2K\~-1/2 ass 
S_. —_— — 
aa = garg (ize) ( Ka )( Kee ) E& Gal Ai (ae) 


The derivative is not a surface term since 8,5 contains tan Q,., in fact it 
can be expressed as 


ds\ 1 qi sin 2Q,, 
(59.) = zor desin 20,2 pene ye See ea} 


qe sin 2Q/, 
ree. _ lo? §=(142K)?—303(1+3K+2K?)+oAK 


The former is a surface and the latter a bulk term. 
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Similarly, for antisymmetric vibrations 


Ay f KO\? (142K 142K\-12 05 4 
Alene ea yf oS ees > 
vt aara(te) (ae )(-Ter) [+ GR), ] oo 


where 


06 A i) de: q? sin 20, 
— aa eer ELD, NM), 3 eA | 
Ce), i m(a) 2 sin 20, 0 A eae (47) 


The surface term of the total density p,=p,4+p,% thus disappears and 


Pan A Le (dee) (orm 2K (gf rin 20, 
Pi oP 2\1LK Ka Kee QP sin 2Q,, 
Using the alternative definition of density (43) involving 6, instead of 5, 
gives 


A K \in 14+2K\-12 qf 1sin2Q, 
P2—= 5-37 2 aT LL |—>- . . . : (49) 
27D? \I+K Ka gq; «sin2Q ,, 
Although the two values of p are formally different, they give the same 
result when integrations over the Q,., Q,. plane are carried out. This is 
due to the fact, that when integrating over one of the small squares of area 


7/4 in fig. 2 the variation of the algebraic terms is negligible compared 
to that of the trigonometric terms. Using (6) 


(pat 1/2)n gin 2Q dQ ger g sec? Qn 
Sp SIO meaty 1+g? tan? Q,, 


where 7 is half an odd or even integer. Thus the integral over the small 
squares is —7?/4. The same result applies for the antisymmetric solution 
(7). Thus both p, and p, have the ‘ average ’ value (taken over the small 
squares) 


(48) 


dQu=—5, - - (50) 


A K \ir 2K +1\-12 142K : 
ih ae laa ban 
p sae (Gan) (: Ke ) (14 Ko ) (01) 
It follows that on the average 
sin 20, rt 
a Tape mee Qian so. 6 (52) 


where n is an integer. In fig. 2 we also exhibit the solutions of (7). The 
two sets of solutions may be shown to be approximately symmetrical 
about the lines given by (52). Now using (52) as the frequency equation 
would lead to a density given by (51). Thus it appears, that the indepen- 
dence of the density on the surface is due to the deviations from the bulk 
distribution, for symmetric and antisymmetric waves, being equal and in 
opposite directions, while those for the waves of type (2) are equal and in 
the same direction. 
The total number of modes is thus 


a2 -L—Pyy?/cos 29}? ay 
% ~ 0 


where @=Q,7+@,?; tan ¢=«, 
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no zero wave number correction being required. Therefore 


A Im r Im : — K ns 
Mplery)= 5 | 7 + a tan “(x Sore gh OE 


where ¢,, =©,,/¢, is the maximum longitudinal wave number. Thus finally 
adding the pontantnion: for the three types, the total number of modes 


M(w,,)=A [fe af a Yt arp +4-I, ee (55) 


Aq 
Substituting the values for K=0, the result agrees with that in B (except 
for the } due to the zero correction) which is obtained from types | and 2 
alone. 
Before calculating the Debye frequency and specific heat of the lattice 
we extend the treatment to include three-dimensional waves and 
recalculate WM. 


§3. THE THREE-DIMENSIONAL CASE 


In the y direction we now take an additional periodic boundary condition 
with periodicity 2L,, which introduces a factor exp (iq, y) in the displace- 
ment, where 

dy nyt Ly§ Ny=.. «. — 38, —2, 15.0, Looe) ade eo) 


Introducing the vector 


q=i+¢,); Q=4q 0, - ei) ey 
where i, j, and k are unit vectors, in the a, y, and z directions respectively, 
it can easily be shown that the two-dimensional motion discussed in the 
last section takes place in the plane defined by q and k. (q replacing 
q, m the various formulae.) However there is an additional type (4) 
consisting of purely transverse bulk waves whose displacement is perpen- 
ducilar to the plane 4, k. (These are omitted in B.) They are again 
symmetrical or antisymmetrical with frequency equations given by 


Q = yn/2 > re= { 


teplacing the parameter Q,, of the family of hyperbola (10) by Q gives 
the number of these for a range dQ of Q as 


1, 3, 5, ... symmetrical, 


2,4,... antisymmetrical. 


> 


pi) am Hater ag _ ee g tities t 59) 


Using this result the densities and total aia eis of modes for the various 
types are : 


(1) "(Os y= ge a ee Oe OE SERE 
> 2K+1 ; 6! Ka” 
i ers RT) [i (ie oe Oe 
V 1+K \1?2 l 1 
h — | imi Gn aT , 
My(w m)= qe am 6 (oe) + a (= = z Js) | : ‘ . (62) 
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where 
74/2, 
1K) |. (1+ 5% Tok tan? 8) 8 (Gato) dhe ee. 9 (63) 
V Ka ae 

Tay fe Lap e A 
ey [te | (qs) an 8 
(4) pee el tae Ve ee, el, 66) 

Gt? Gar A 
M pm) =v | ts Gigs Ts SelM ee Ore teke ee ECT) 
Thus summing over the four types the total number of excited modes is 

9 Int Us Ger. 

Wm) =V [i a aos tad, (Y,? +4 De 1) | : : . (68) 


$4. THE Sprcrerc Hear awn ENERGY SURFACE TERMS 


The Debye frequency can now be determined by equating the number of 
degrees of freedom (3N V) of the lattice points to the number of excited 
modes M(W,,,), hence 


BN = (2G no Imi) OT + Imel Ve +h—le)f4nk, . . (69) 
Thus when L,>00; Gin Gnt—> Imo dnt where 
Doe Peg Pee Br Ne 1 |, P8070) 
giving the Debye frequency and temperature (4p) for a large specimen, 
w= b6r2NC or Op =hw,,|k=he(6r2=N)V3/k . 2. (71) 
where the ‘ mean velocity ’ c is defined by 
Bie WGP ieee ts ts ew (2) 
Substituting hag for q,,, in the surface term of (69) gives 
ponies ON eT. 8. 2, (73) 


where 
m(Yot3—Ls) 

218m?) F[2+ [K/(14 2K) pepe 

To calculate the intrinsic energy H of the vibrations we take (68) as a 
general formula for excited modes whose frequency lies between 0 and w. 
The number of modes whose frequency lies in the range w, w-+dw is 
(dM /dw)dw. It is however, simpler to consider ¢ as a function of q,(=«/¢;). 
Thus 


Wik 


Unt aM hdc; = 
“= — d he = —__—__..._.. _ (4 
j, <(d0) dq, bE uae exp (hiq¢,/kT)—1 


Evaluating this integral with the usual new variable 
wahge|bl; e=hq,cfkT—Oy/T . . . . (75) 
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and substituting from (73) for the upper limit 


bes a ea ae fof 


9 et — ic L, 


The first term represents Debye’s result while the second differs from that 


given in B only in the constant A(4). 
We have evaluated the various constants for three values of K. The 


results are collected in table 1, where the corresponding Poisson’s ratio (a) 
is also shown. 


K= 


Table 1 
K y I; Xo Jig: (Yo?+3—Ls) A 
0 4 1-35 0-2956 1-095 0-245 0-044 
1 4 0-836 0-4641 1-18 0-84 0-14 
oe) 0 0-588 0-6181 1-31 1-22 0-19 


The single value of the constant A(K) given in B is %0-24 which is 
about five times as large as the corresponding correct value A(0)0-044. 
However if we use elastic constants corresponding to A=1, the results in 
B will only have to be scaled down by a factor of about 0-6. 

From (76) we can obtain the surface free energy Fg of the lattice 
vibrations 


P,=E—T [7 7 (an) ¢7=— 9A(K)kOpN2°Z (x) . . (77) 
1 (* udu 
where Zit) = ea ‘ (78) 


This result, with the constant as given in B, is used by Huang and Wyllie 
(1949) in their paper on the surface energy of a metal for the temperature 
dependent energy contribution of the lattice. 

Finally the specific heat can be written in the form 


Op=Of1+A(K(myNL,] 2 . 2°. . (79) 
where Ry eed Cae a 
x 
all eae 


and C’, is the Debye specific heat. An extensive discussion of this result 
is given in B where curves are drawn to show the variation of C,, with 
temperature. These must be corrected by the scaling factor A(K)/0-24. 


. CONCLUDING REMARKS 


It is interesting to ae that the reason why our results are in 
formal agreement with those of previous calculations is that the number 
of modes M(W,,) has a surface term proportional to mt. (OT Fmi)-28 for 
example in (68). It seems unlikely that more claborate ‘boundary 
conditions would alter this situation, Their only effect would be to 
modify the constant A(K). 
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If we consider a cubic lattice then q,,~4/a where a is the lattice spacing 
(taking c,=c,). Thus the condition that the total energy of the specimen 
can be split into a bulk and a surface term, i.e. g,,/,> 1 is equivalent to one 
of the conditions for replacing the lattice by a continuous medium, viz. 
a<l,. Thus if the Debye theory is applicable to a large specimen of a 
material then the above calculation will give surface terms for smaller 
specimens to the same degree of applicability. 

The only experimental evidence having a bearing on this work is the 
difference in the specific heat of activated charcoal and graphite observed 
by Simion and Swain (1935). This is attributed to the greater specific 
surface of the former. In B the specific surface of the activated charcoal is 
estimated to be about 500-800 m?/g. This area will be slightly increased 
when our values for A(X) are taken which does not affect the interpre- 
tation of the result. 
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APPENDIX 


We briefly show how far the extension from two to three dimensions 
carried out in B is not quite correct. Since the authors use a different 
representation of the modes of vibration we first repeat our own calculation 
in their representation for the special case K=0 treated by them. 

In B, representative points are considered in the q,, 7,, Y, Space instead 
of the Q;,,, Q’;, plane for the case K=0 when q';,=q,. This would be 
inconvenient for other values of K. 


Thus ; 
nil 0, = al and.g,—,,m 21,0), . . (80) 
where tan 6=¢(«’); 0<d<n/2 giving the density 
V 5’(a’) 
== Sep) (tie 2 Vali Meare teal 
named ae (1 OL, ) a € (81) 


considering symmetrical and antisymmetrical waves together. 
Introducing polar co-ordinates with 0q,, as axis (q,, 0, %) 


tan 0=q,/¢,; tanv=lWo'=ltang. . . . . (82) 
Vi rl 2% me d’(tan¢)\ : 
M (0m) = el, J, I, (1 = q,L, cos d q¢ COS p dd do dq, ea. (83) 
mt _— Int RF ee aeae (84 
=v | ed ars kote: (84) 


which agrees with (62) (K=0) after adding the Z, correction, 
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In B, the extension to the three-dimensional case is carried out by adding 

a factor exp(iq,y) to the displacement function for the x, z plane. Thus 

x’ is not taken as the tangent of the angle between q,, and ¢(=q’,,) but 

of the angle between q,,andq,. The authors then take polar co-ordinates 
q'4', d' with 0q,, as axis 

cos 0’=q,/q,, tan ¢’=¢1/4e- ce St es Me 


This gives the wrong result 
pena V m/ 2, 2n Imt 3’ (tan ¢’) : , 
IV = — = 2 gin A’ dé’ dd'd 
Vy (w,,) Zz i. iF if (1 mS vr) q; sin 0’ dé’ dd’ dq, 


Mylon)=V Yt Imi (0) (86) 
M m 6772 2n?L, 2 . . . . . 


which agrees with their eqn. (29) 


po 
Mi(r) = TP lbs 7 


when account is taken of their notation 
r= (mil )(v/n) 18 and. n=L,/L,. 
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SUMMARY 
By the accurate control of the energy of a 20 Mev betatron it has 
been shown that the activation-energy curve consists of a series of 


small plateaux indicating the existence of narrow resonances for the 
photo-neutron reaction in copper. 


=e 


$1. INTRODUCTION 


In recent years numerous publications, e.g. McElhinney ef al. (1949), 
Johns et al. (1950) have appeared describing work carried out on 
photo-nuclear reactions. The usual technique involved in these 
experiments consists of exposing a suitable foil to the high energy 
bremsstrahlung from a betatron or other suitable electron accelerator. 
_ Many of the nuclei produced in these reactions are radioactive. This 
activity serves as a measure of the photo-neutron yield and by varying 
the peak energy of the incident photons an activation curve can be 
deduced. 

By making certain assumptions about the tip of the bremsstrahlung 
spectrum a rough estimation of the (y,n) cross section can be made. 
Previous authors, however, seem to have been content to vary the 
peak energy of the accelerated electrons by intervals of the order of 
0:25 Mev. This only indicates the general shape of the excitation curve 
(see fig. 1) and does not bring out any fine structure which may exist. 

This report describes an experiment in which the peak energy 
of a 20 mev betatron has been varied by 0-02 Mev intervals. The 
®8Cu(y, n)*Cu reaction has been taken as the subject of investigation 
firstly because ®Cu has a reasonable half life (approximately 10 minutes) 
and secondly as this reaction has been investigated so many times 
previously it has become to be almost a standard reference cross section. 


§2. ConTRoL OF THE BETATRON ENERGY 


The kinetic energy of the accelerated electron in a betatron is 
determined by the magnetic field at the orbit and the radius of the 
orbit. The peak x-ray energy is then dependent on the excitation 
current which passes through the magnet coils, that is provided the 
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expansion of the electron orbit takes place at maximum magnetic field. 
The magnet current can be controlled manually within a few per cent 
by means of a buck-boost circuit. A finer control of the energy can be 
further achieved by expanding the orbit at a fixed value of magnetic 
field. This has been accomplished by using a current transformer 
connected in series with the magnet. This transformer provides a 
voltage, proportional to the magnet current, which is fed to a multiar 
circuit which produces a square pulse having a front of extremely short 
duration. The multiar was developed by Williams and Moody (1946) 


Fig. 1 


Relative activity (counts/roentgen). 


Peak energy of betatron y-rays (Mev). 


A typical activity curve for the °Cu(y, n)®Cu reaction. 


for use as an accurate range marking device in radar. The circuit is 
simply dependent on the conduction of a diode which only occurs when the 
voltage of the controlling sine wave has achieved a predetermined value 
fixed by a stabilized voltage, fig. 2. The voltage can be precisely 
adjusted by the variation of an accurate potentiometer (nominally 100 kQ) 
and by this means it can be seen that the multiar pulse can be made to 
move up or down the sine wave. Furthermore, provided the current 
transformer has a zero phase shift the multiar triggering position is 
independent of the magnet excitation. The square multiar pulse is fed 
into a pulse shaping circuit which produces a standard trigger voltage 
This in turn fires the hydrogen thyratron on the expander modulator. 
and also operates a megavoltmeter similar to the one described by 
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Westendorp (1946). The megavoltmeter integrates with time the voltage 
induced in a loop around the betatron pole pieces, between the limits 
imposed by the times of injection and expansion. This integration 
produces a voltage proportional to the product of the field and the orbit 
radius, i.e. the peak x-ray energy, giving an indication of the electron 
energy with an accuracy of approximately 0-20 mev. Additional 
checks on this energy were provided by the monitoring C.R.O. and also 


Fig. 2 
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A block diagram of the energy control circuit. 


by measuring the d.c. stabilized voltage fed to the multiar circuit. The 
methods of energy determination although helpful in carrying out the 
experiment cannot be expected to give an accuracy of the energy value 
better than +1°% and have not been relied upon in the present work. 


§3. Discussion OF THE ACCURACY OF THE METHOD 


As already stated the only two factors governing the control of the 
firing of the multiar are : 

1. The time of firing of the diode in the multiar circuit which is 
controlled by the d.c. reference voltage. 

2. The phase shift in the current transformer. 

An attempt to control the first factor was made by using a stabilized 
power pack which employs the usual series tetrode method providing a 
stability of 1 part in 5000 at 500 volts output for a wide range of input 


an? 
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mains voltage. A further improvement was obtained by the use of a 
constant voltage transformer ; all the electronic circuits obtained their 
50 cycles supplies from this transformer. 

The phase-shift of the current transformer was made small by careful 
design, but the first attempt turned out to be somewhat of a failure. 
The phase shift-excitation characteristic was measured carefully and 
indicated a phase shift of two degrees at 60 amp falling gradually by a 
total of 30 minutes at 120 amp. This phase variation of half a minute 
per amp was good enough for relative energy measurements and many 
of the preliminary energy checks were carried out using this transformer 
and keeping the magnet excitation current constant to within a few amp 
by manual control. The second and more successful transformer, now 
in use, was constructed of mu-metal laminations and gave a total phase 
shift of 15 minutes over the same current range and with a correspondingly 
smaller relative change. 

Another factor which may cause variation in the electron energy is 
the delay between the firing of the diode in the multiar and the triggering 
of the expander. This time delay has been accurately measured and 
accounts for a: possible increase of 0-01 Mev in the final energy. This 
error does not affect relative measurements in energy but only absolute 
values which can be suitably corrected. Summing up, it has been 
estimated that the relative energy values can be held to within 
+0-02 Mev. 

§4. EXPERIMENTAL DETAILS AND RESULTS 

The layout of the experiment is shown in fig. 3. The foils made of 
1 in. square x 0-060 in. thick spectroscopically pure copper were exposed 
at a distance of 30 cm from the x-ray targets. The foil itself was fastened 
by adhesive tape on to a bakelite jig thus enabling similar foils to be 
placed in an identical position. Some previous workers have tended to 
place the foils as close to the target as possible. This reduces the 
exposure required for a given count, but it places a severe restriction on 
the position of the foils owing to the small size of the x-ray beam, and 
any variations in foil position would be likely to cause errors. 

The foils were exposed for several minutes and the integrated x-ray 
intensity was measured by means of two 25R Victrometer thimbles 
placed in a 6-in, diameter polystyrene cylinder placed at 75 em from 
the target. Two thimbles were used in an attempt to avoid any errors 
in measurement and if the ratio of the readings of the two thimbles was 
not within 2% of the mean value of the x-ray dose the exposure was 
repeated. Several runs were carried out at each energy point in order 
to obtain consistent results. The betatron energy was varied by steps 
of 0-02 Mev. 

The activity of the exposed foils was measured by means of an 
end-window Geiger counter and the usual electronic scalar. The 
exposure times were controlled so that the numbers of counts were 
greater than several thousand per ten minute period, thus reducing the 
probable error to less than 1%. The counter equipment was housed 
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in an area well away from the betatron enclosure, enabling counting 
to be done while other exposures were being carried out. This procedure 
was necessary in order to reduce the time of the experiment as several 


Fig. 3 
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A schematic diagram of experimental arrangement showing the location of the 
copper foil and the Victrometer thimbles. 


Fig. 4 


Relative activity 


Maximum y-ray energy (Mev). 
The activation curve for the energy range of 13-5 to 17-5 Mey. 


thousand separate exposures had to be made. The plot of foil activity 
per Roentgen is shown in fig. 4. It can be seen that the excitation curve 
is of the same general shape as fig. 1 but it actually consists of a series 
of sharp rises and horizontal plateaux indicating the presence of multiple 
resonances.* The width of the observed plateaux decreases from 
approximately 0-4 Mev at an excitation of 14 Mev to about 0-08 at 


* A similar effect had been reported recently by Goward and Wilkins (1952) 
in their studies of the C(y, 3a) reaction in photographic emulsions. The 
resonances are much broader as would be expected for a lighter element. 
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16 mev. ‘There are most likely additional resonances in the higher 
energy region which cannot be resolved by this method. Above 16 Mev 
the data hardly warrants drawing a curve but there is some indication 
that other small steps exist. At the lower end of the curve the statistical 
inaccuracies increase because of the lower cross section and the 
corresponding lower betatron yield (the Roentgen output drops by a 
factor of about 10 when the peak energy is reduced from 20 to 10 Mev). 
For these reasons the steps become indistinct near to the threshold of 
the reaction. The energy scale in fig. 4 was deduced by the calibration 
of the multiar potentiometer against the ®Cu(y,n)®Cu threshold at 
10-9 Mev. It must be borne in mind that this threshold is not known to 
better than -.0-2 Mev and consequently the energy scale may be in error 
by as much as + 0-25 Mey. 


$5, THEORETICAL CONSIDERATIONS AND DISCUSSION OF THE RESULTS 

In the experiment described in the previous sections the activity of 
the copper sample is given as a function of the peak energy of the 
bremsstrahlung, or the energy of the electrons accelerated in the 
betatron. In order to interpret fully the shape of these curves it is 
obviously necessary to know the number of quanta produced in the 
forward direction. This energy distribution has been deduced by 
Schiff (1946) from the theoretical work of Heitler (1944). The experi- 
mental spectrum has been determined by several workers, see Koch and 
Carter (1950) and in particular by Phillips (1952) for the betatron used 
in this experiment. The results in the main correspond to the theoretical 
results but they suffer from poor statistics at the higher energy end and 
do not indicate the shape of the curve in that region. 

No work seems to have been reported on measurement for the 
higher energy end of a spectrum in the region of 20 Mev, but Miller 
and Waldman* (1949) describe the result of the investigation of 
bremsstrahlung by means of activity measurements on the nuclear 
isomerism of Indium for the energy range 1-0 to 2-6 Mev. Furthermore 
the theoretical results are only approximate at the high energy and 
since the calculations are based on Born’s approximation, which requires 
that 27(Ze?/hv) <1 and 27(Ze?/hvy) <1, where Z is the nuclear charge of 
the x-ray target, e is the electronic charge, vy and v are the velocities of 
the electron before and after the creation of the radiation. It can be 
seen that even for high energy electrons and elements of high atomic 
number this condition does not hold, e.g. for platinum Ze2/fic=0-54. 
(A theory of high-energy bremsstrahlung has been produced by Bess 
(1950) which avoids any approximations. Unfortunately only one point 
on the curve is given and the numerical calculations involved are rather 
intractable.) These approximate results show that the number of 
quanta falls to zero at the high energy limit. The more exact theory of 
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Sommerfeld and Maue (1935) shows that for non-relativistic energies 
this is not the case and the number of quanta is finite. Guth (1941) has 
argued quantitatively that this exact non-relativistic formula may be 
obtained from the results given by Born’s approximation by multiplying 
them by a function f(Z, v, vy). On these grounds it seems reasonable 
to use a similar function to convert Heitler’s results into a formula 
which will provide a good approximation to the actual values. For 
interpretation of the experimental results it is advisable to introduce 
the conception of an isochromat. This is defined as the intensity of a 
narrow band of x-rays of a given energy hyp as a function of the peak 
X-ray energy (electron energy). The isochromat for the Heitler theory 
turns out to be a curve steadily increasing with energy as shown in 
fig. 5 (a), whilst the corrected energy distribution has the form of a 
unit function as in fig. 5 (6), i.e. the isochromat is independent of energy. 


Fig. 5 


INTENSITY INTENSITY 


PEAK ®’T = y-RAY PEAK “Vr y-RAY 
ENERGY ENERGY 
(a) (d) 


(a) The isochromat deduced from the Heitler calculation using Born’s 
approximation. (b) The isochromat deduced from the more exact non- 
relativistic Sommerfeld theory of bremsstrahlung. 


From the shape of these isochromats it is easy to show that for a 
steadily increasing isochromat one would expect an activity function 
consisting of a series of sloping lines with discontinuities at the various 
energy levels, see fig. 6 (a). For a step like isochromat an excitation 
curve having a series of sharp rises and plateaux would be expected, 
fig. 6 (b). A slight slope of the front edge would occur owing to the energy 
spread of the accelerated electrons. The third curve, fig. 6 (¢), indicates 
the activity function for thick target radiation. This can be understood 
by assuming that the thick target isochromat is the sum of a series of 
thin target isochromats thus producing a sloping activity curve with 
discontinuities at the different energy levels. 

Summing up it may be said that this work indicates : 

1. That narrow resonances exist in the ®Cu(y, n)®Cu reaction and 
that the cross section therefore should be deduced on the basis of energy 


540) K. Phillips on the Existence of M ultiple 


levels in the compound nucleus. This does not eliminate the possibility 
of a photo-electric effect and it may be concluded that the resulting 
cross section is due to a series of narrow resonances superimposed on a 


continuous background. 
2. The high energy bremsstrahlung spectrum 1s of the thin target 


type, thus confirming previous spectrum measurements. 


Fig. 6 


ACTIVITY 
ACTIVITY 


v-RAY ENERGY PEAK -y-RAY ENERGY 
(a) (8) 


ACTIVITY 


PEAK ‘Y-RAY ENERGY 
(c) 
(a) The activity curve deduced from a steadily increasing isochromat. (6) The 
activity curve for the case of a step isochromat. Both (a) and (6) are for 


thin target radiation. (c) The activity curve for a thick target radiation 
for small thin target steps. 


3. The theoretical bremsstrahlung curve of Heitler’s is unsatisfactory 
at the high energy limit owing to the use of the Born approximation, 
and furthermore the isochromat has the form of a step function. It is 
worth while noting from fig. 4 that the front of the step is not a sharp 
discontinuity at each level. As already mentioned this might be caused 
to a certain extent by the energy spread of the accelerated electrons. 
From the curve the average energy variation over the rapid rise is 
approximately +-0-04 Mev—slightly larger than the previously estimated 
energy spread of +-0-02 Mey. 
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§1. INTRODUCTION 


Iv is rapidly becoming clear that there exists a very large measure of 
charge-independence in the interactions between nucleons. This reali- 
zation has given rise to a revival of interest in the isotopic spin formalism 
with its provision of another significant parameter for the description of 
nuclear states. Adair (1952) and Kroll and Foldy (1952) have discussed 
the effect of conservation of isotopic spin and charge parity in heavy 
particle reactions ; Trainor (1952) has pointed out the existence of isotopic 
spin rules in the emission of electric dipole radiation and Radicati (1952) 
has shown that, for any multiple whether electric or magnetic, the rigid 
selection rule 47’=0, +1 obtains, with the limitation that for electric dipole 
transitions in self-conjugate nuclei the selection rule is 47=+1. Evenif 
the specifically nuclear interaction between nucleons is charge-independent, 
complete charge-independence cannot obtain because of the existence of 
Coulomb forces; we must therefore expect that the isotopic spin is not a 
good quantum number and that, for example, there will be a mixing of a 
small amount of a wave-function belonging to a state of 7’=1 ina state that 
is principally of 7’=0. This will effect a partial relaxation of the isotopic 
spin selection rules; if they are obeyed as well as may be expected, 
having regard to the Coulomb perturbation, the hypothesis of charge- 
independence of the specifically nuclear forces is strengthened, if they 
appear to be poor rules charge-independence must be only approximate. 
Alternatively, if we can obtain some independent assurance as to the 
validity of charge-independence such as may derive from nucleon and 
meson scattering experiments, violations of the isotopic spin selection 
rules will afford quantitative information about the properties of groups 
of nuclear levels which are not directly investigated. Radicati (1953) 
has investigated theoretically the mixing of isotopic spin states by 
Coulomb forces in the ground states of certain nuclei (with two and four 
nucleons outside a closed shell), and finds that, in his examples, we may 
expect to find mixed with the principal wave-function appropriate to 
T=0 a few tenths of a per cent to a few per cent (in amplitude) of a wave- 
function appropriate to 7'=1. This means that, under complete specifically 
nuclear charge-independence, we should find that transitions that are 
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wllowed by ordinary spin and parity rules but forbidden by isotopic spin 
rules occur about 10? to 104 times less readily than those allowed by isotopic 
spin rules as well.* 

This paper reports experiments that place rather crude upper and 
lower limits to the validity of the isotopic spin rules; these limits are 
In agreement with the above-mentioned work of Radicati. 


Fig. 1 
7116 i 
6-913 tt 
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Low-lying levels of 1°O. The full arrows indicate known transitions, the dashed 
arrows those sought in the present investigation. 


§2. THE Low-Lyine States or 16O 


The present investigation concerns itself with the rule that a transition 
with A7=0 is forbidden for electric dipole radiation in a self-conjugate 
nucleus; 16O appears to provide examples both of the operation and 
violation of this rule and so enables us to set rough limits to its validity. 

We consider the lowest states of 1°80; these are shown with their 
relevant characteristics in fig. 1. The energies are taken from Ajzenberg 
and Lauritsen (1952); the spin and parity of the first excited state 


* The mixing of isotopic spin states will depend on the proximity of the 
‘contaminating’ states to that in question and so, presumably, be least for 
the ground states ; the rules, therefore, should become apparently less strict 
for transitions involving highly-excited states. 
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are due to work by Rasmussen, Hornyak, Lauritsen and Lauritsen (1950), 
Devons and Lindsey (1949), and Devons, Hereward and Lindsey (1949) ; 
those of the (3—) state to Arnold (1950) and Barnes, French and Devons 
(1950); those of the remaining two excited states to Seed and French 
(1952), Sanders (1952) and Fagg and Hanna (1952). The known modes 
of de-excitation are shown by the full arrows ; the hypothetical modes of 
interest in the present investigation are shown by the dashed arrows. 
Now 160 is a very strongly-bound nucleus and the ground states of its 
neighbours 1*N and !8F (the latter unstable to break-up into 1°O and 
a proton) occur, after allowance for the Coulomb energies and n—p mass 
difference, 12:8 Mev above that of 16O (Lauritsen 1952, Ajzenberg and 
Lauritsen 1952). It is therefore certain that all five states of '*O shown 
in fig. | have 7'=0; if we think in terms of charge-symmetry alone, the 
charge parity of all these states is almost certainly even. 

It is apparent that y,, which must be electric dipole, violates the isotopic 
spin rule.* If this rule were strict the (1—) level would de-excite by 
means of the cascade starting with the electric quadrupole y,; we have 
sought ys, as will be described later, but have not found it. The (2+) 
level emits the electric quadrupole y, although we should ‘expect the 
cascade starting with the electric dipole y, to compete favourably with it ; 
we have sought y, without success and its absence may be an illustration 
of the isotopic spin rule in action. The quantitative implications of 
these observations will be discussed after the account of the measurements. 


§3. THE POPULATION OF THE INITIAL STATES 


In order to search for y, and y, we must be able to prepare !6O in the 
(l—) and (2+) states. The most convenient method is to bombard 
fluorine with protons when the reaction 1F(p«)!®O takes place. This 
reaction is strongly resonant below I Mey of proton energy and by 
choosing the appropriate resonance we may vary the population of the 
states shown in fig. 1. Unfortunately, all resonances below 1 Mev that 
populate the (1—) and (2+) states favour the (3—) state, which merely 
provides the useless y,; as background. We have used the resonances 
at 874 kev and 935 kev proton energy ; the former populates the levels 
(1—), (2+) and (3--) in the ratios 7:19:74 while the latter gives the 
ratios 22: 3:75. (Chao, Tollestrup, Fowler and Lauritsen 1950, Seed 
and French 1952.) 

$4. Tor MEASUREMENTS 


We have examined the gamma-rays directly using Nal(T1) crystals ; 
the whole gamma-ray spectrum was examined to begin with to make sure 
that no medium-energy gamma-rays comparable in intensity with y,, ys 
and y; existed. We then examined the low-energy region of the spectrum 
with some care using a kick-sorter of the type described by Hutchinson 
and Scarrott (1951). Several examinations were made with good statistics, 
See Se eS ee 

* This has also been pointed out, in conversation, by Professor W. E. Burcham. 
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achieving the order of 100,000 pulses per channel in each one. Each 
examination was split into many runs, between which irradiations with 
'87Cs gamma-rays enabled a check to be kept on the overall stability of 
the system. We have found it best to display graphically not the actual 
counts per channel, since the very small statistical errors do not permit 


Fig. 2 
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Differential distribution of NalI(TI) crystal pulses in successive channels at 
proton bombarding energy of 874 kev (thin target). The arrows 
indicate the expected positions of the various low energy gamma-rays 
(see fig. 1). The full line shows the expected differential distribution due 
to a branching ratio of 5x 10-*: | for y, relative to yz. The error limits 
represent standard deviations. 


of any easy demonstration of the sensitivity of the method, but rather the 
differences between actual recorded counts in successive channels ; this 
gives a differential distribution. Figure 2 shows such a differential 
distribution for a single examination of the spectrum at the 874 kev 
resonance, using a thin target and the detector at 90° to the proton 
beam.* The arrows indicate the positions at which the various low-energy 


SO Ie AS Ee Se ee ee eee 

*It is very unlikely that strong angular distributions obtain at either 
resonance ; those for the high-energy gamma-rays are known. to be feeble 
(Devons and Hine 1949, Sanders 1952, Chao 1950). 
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vamma-rays of fig. | are expected to occur and it is obvious that nove is 
present in excess of the statistical fluctuations. As an illustration of the 
sensitivity of the method we have indicated by the full line the differential 
distribution to be expected if the branching ratio of y, relative to ys 
were 5X 10-3: 1. Other examinations were consistent with this conclusion 
that ys is at least 200 times stronger than y4. The branching from the 
(1—) level is best sought by bombardment at 935 kev ; the results were 
similar to those of fig. 2 and suggested that the branching ratio of y, 
relative to y, is less than 8x 10-3: 1. y, which might lead from the (1—) 
to the pair-emitting (0-++) level is not found. The experimental branching 
ratio is less than 10-2: 1; the ratio expected from the ratio of the energies 
of y, and y, is about 3x 10-3: 1 so the absence of y, is not noteworthy. 
The present method is not the best available for seeking y,, and we hope 
to make it the object of a separate investigation. 


$5. DIscUSSION 

In order to discuss the significance of the absence of y, and y, we must 
have some idea of the relevant branching ratios to be expected “in the 
absence of isotopic spin rules’. It is not possible to obtain more than a 
very rough orientation on this point and this we derive from the predictions 
of the single-particle model (Weisskopf 1951, Blatt and Weisskopf 1952). 
This model is clearly not directly applicable to the present nucleus, but 
its predictions will at any rate provide an order-of-magnitude estimate. 
It appears that the predictions of this model are rather reliable for 
electric dipole transitions in light nuclei (Wilkinson 1953). 

According to the single-particle model y, should be 2:4 10% times 
more probable than y, ; we find it to be at least 120 times more probable 
so the isotopic spin rule has not slowed up the transition by as much as a 
factor 2 10°, i.e. if we take the ground state of 1*O0 to be pure 7'—0, 
the relative amplitude of the 7’=1 state admixed to the predominantly 
T=0 state of 7-12 Mev is greater than 2x 10-3. 

According to the single-particle model y; should be 6 times less probable 
than y,; we find it to be at least 200 times more probable so the isotopic 
spin rule has slowed up y, by at least a factor 1200, ie. the relative 
amplitude of the 7'=1 states admixed to the predominantly 7'=0 states 
at 6-91 Mev and 6:14 Mev is less than 3 10- (assuming. equal ‘ con- 
taminations ’ of the two states). 

If, then, we may give our results their most naive interpretation it 
appears that these isotopic spin or charge parity states are pure to within 
somewhere between 0-2°, and 3% in amplitude :* this agrees very well with 


* The principal factors likely to disturb seriously our present estimates of 
the expected ratios ‘ without isotopic spin rules’ are correlations which would 
inhibit electric dipole transitions and the possible relevance of the alpha- 
particle model (Dennison 1940, Seed and French 1952) which would do the same. 
It is unlikely that the former effect is serious for so light a nucleus with an 
excitation of 7 Mev (see Wilkinson (1953)); about the latter little can be said 
at the present time. If either of these effects is of importance both limits must 
be increased ; they are probably correct to within a factor of about five, 
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Radicati’s estimate of a few tenths of a per cent to a few per cent as the 
contamination in amplitude due to Coulomb forces, and so may be taken as 
some strengthening of the hypothesis of the charge-independence (or 


charge-symmetry) of the specifically nuclear component of the nucleon— 
nucleon interaction. 
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ABSTRACT 


It-is shown that an ideal Bose-Einstein gas confined to a layer of any 
finite thickness behaves like a ‘ two-dimensional gas ° in that the tempera- 
ture of the B.E. condensation is size dependent and tends to zero as the 
area of the layer tends to infinity. Such an effect is not observed in real 
helium films. If it is assumed that liquid He must be treated as if made 
up of independent B.E. assemblies, each smaller than a fixed length 
Ly (~7xX10-&em), the observed variation, with thickness, of the 
temperature of onset of superfluidity, in unsaturated films, may be closely 
predicted. The zero-point motion in each restricted assembly is then 
close to the critical transfer velocity of superflow, and L, is of the same 
order of magnitude as the thickness of the saturated film. This assump- 
tion brings the ideal B.K. gas model nearer to the real structure of liquid 
He, but the interpretation of L,) as a ‘ scattering length ° does not give 
agreement with estimates from kinetic theory. 


§ 1. INTRODUCTION 


Ir is well known that one interpretation of the lambda-point transition 
in liquid helium is that it is equivalent to the condensation process which 
can be shown to occur, theoretically, in an ideal Bose—Kinstein gas. 
The details of the process are not well predicted, but the absence of any 
similar phenomena in *He, which obeys Fermi—Dirac statistics, seems 
very strong evidence that this hypothesis is fundamentally correct (see 
review articles by Dingle 1952, Daunt 1952). The mathematical argument 
for the condensation has been exhaustively studied by numerous authors, 
and there now seems little doubt that in an infinite assembly of B.E. 
particles there is a discontinuity in certain thermodynamic properties 
at a definite temperature. It is usually assumed that a macroscopic 
sample of liquid helium may be treated as if it were part of such an 
infinite assembly, and that the appearance of superfluidity corresponds 
to the condensation process. 

There are serious difficulties, to which we shall return later, in this 
treatment of a liquid as if it were a perfect gas. For the moment we 
observe that the phenomena of superfluidity are not confined to the 
bulk liquid, but have been observed, at sufficiently low temperatures, 
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in films only a few atomic layers in thickness (Bowers, Brewer and 
Mendelssohn 1951, Long and Meyer 1953). If the assumption is valid 
that the onset of siportuidity corresponds to a B.E. condensation, this 
should be true as much for a thin film as for the bulk liquid. It is, 
therefore, of some interest to study the mathematical properties of an 
ideal B.E. gas, not, as in the usual theory, in an infinite cubical box, but 
confined to a layer of finite, even microscopic, thickness, and to see 
whether these properties agree in any way with the observed properties 
of thin films of helium. 

If we treat a film as a two-dimensional gas, which would, at first sight, 
seem a reasonable approximation, we find that there is no B.E. conden- 
sation at all, according to the standard argument, in an infinite square. 
For a finite assembly of this type, Osborne (1949) has shown that there 
exists an ‘ accumulation ’ temperature, akin to a condensation tempera- 
ture, and varying as 1/In VN, where N is the number of atoms in the 
system. 

The problem of the ideal B.E. gas in a layer of finite thickness has also 
been discussed by Osborne (1949), who has argued that the usual dem- 
onstration of the B.E. condensation breaks down only when the layer 
has become thinner than a critical distance, g, which is given, approxi- 
mately, by the condition that the mean thermal energy, k7', of each 
atom shall be of the same order of magnitude as the zero-point energy 
p*/2M, where p is the momentum conjugate to q, satisfying the uncertainty 
relation pq~h. It would seem evident from first principles that for 
films thinner than qg the higher energy levels, corresponding to momenta 
2p, 3p, etc., normal to the plane of the film, would not be excited, and 
the system behaves as if it were two-dimensional. 

From this argument, one would evidently conclude that if one took 
a film of large, but finite thickness (e.g. 1 cm) and infinite extent, it would 
show an ordinary B.E. condensation, but that as one squeezed the faces 
together the transition temperature would be gradually reduced, and 
would tend to zero when the critical thickness was reached. Such a 
lowering of the onset of superfluidity in thin films is actually observed ; 
the original object of this research was to calculate it explicitly, for it 
appeared that the results were implicit in Osborne’s work, and merely 
required analysis. 

However, it turned out that the difficulty of estimating the errors in 
Osborne’s replacement of sums by integrals led to the adoption in the 
present work of a different technique, by means of which the following 
paradoxical result has been obtained: the condensation temperature of 
a B.E. gas confined to a layer of fixed thickness, however large, tends 
to zero as 1/In NV, when N, the total number of particles, tends to infinity 
(the density being kept constant). In other words, a film of finite 
thickness always behaves as if it were.‘ two-dimensional’. It is impossible 
to assume that a real film is just a small part of an indefinitely large 
assembly; the size of the assembly affects the transition temperature. 
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The obvious way of escaping from this difficulty is to assume that there 
exists a natural upper limit to the size of the assembly, and that macro- 
scopic quantities of liquid helium are to be thought of as made up of large 
numbers of such maximal assemblies, each independent of its neighbours. 
There is some theoretical justification for this device, and it will be shown 
that we can obtain a good estimate of the variation of the lambda point 
with film thickness, without changing the usual results for the properties 
of the bulk fluid. This natural limiting size is a cube of side about 
710-8 cm and it is shown that a fundamental length of this order of 
magnitude plays an important role in many distinct phenomena observed 
in} He Li: 


§2. BosE—-EINSTEIN GAS IN RECTANGULAR Box 


As usual, we take a system of B.E. particles, each of mass M, whose 
wave functions are confined to a rectangular box of dimensions L, L, D, 
where the thickness D<L. The energy levels of the system are then 
given by 


Bin =(P-+-m?)ckT+n2BkT, 2 2... (2) 


where «=h?/SML?kT and B=h?/8MD*kT and 1, m, n, are integers. 
Following Fowler and Jones (1938), we use the method of the Grand 
Partition Function, which tells us that the total number, N, of particles 


is given by 
io) 
N= : a exp (4-9 —JE inn! kT), i ee fide i Tag IK (2) 
j,l,m,n=1 

where kT is the chemical potential of the gas particles. Thus, from 

(1) and (2), we have 
N= 2 expj{u—(?+m?)a—n?B}. Ssh ieee 

glnn 

For a fixed value of j, the remaining terms represent a product of 


Theta Functions, to which we may apply Jacobi’s transformation (Whit- 
taker and Watson 1952, p. 475). This may be written 


<9) foe) 
2 . 2 92 
eed (—ws°)=$4/(n/x){1+2 2 exp (—ms?/z)}-2. . . . (4) 
= &=1 

Now, if «<7/2, the error in taking }{\/(7z/x)—1} in place of the right- 
hand side of (4) is at most of the order of one part in 500, Again, if 
xt >*1/2, the error in taking exp (—2) alone, in place of the sum, is of the 
same magnitude. Thus, it is possible, to a sufficient degree of accuracy, 
to take these simple expressions alone to represent the Theta Function 
for appropriate values of «. 

l'o use this approximation in the evaluation of (3), we must divide 
the summation over j into three regions : 


(a) 1<j <n/28, (b) 7/2B<9 <m/2a, (c) m/2a<4. 
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We here assume that 7/2« and 7/28 are integers; this can easily be 
arranged, since « and 8 are small, except for very thin films. No difficulty 
arises if this condition is not fulfilled; we merely taken the nearest 
integers instead. 


We may now write (3) in the form 


] 7/28 


Ny 2 exp (Hf) mlia)—1}°v/(mliB)—1} 
jJ= 
1 m/e ‘ ; oe) : 
+g % exp titu—A)HV(mija)—1P+ 2 exp {7j(u—2u—)}. 
j=n/2B +1 j=n/2a+1 


(5) 
It is evident that each of the terms arising from the sums over regions 


(a) and (6) is bounded, unless ~ > 0, or «, or B vanishes. But the last 
term has a singularity when 


fore h ii errs 0 Malady eases &. ches (6) 


where the sum behaves like 1/e. To obtain the thermodynamic properties 
ofthe assembly, we must solve (5) for w, as a function of V, etc. Evidently, 
when JN is small, w is negative. As N increases, so does pu, until it 
approaches the value 24+ 8. But here any large value of N can be 
fitted by taking « sufficiently small. On the other hand, the first two 
_ sums are monotonically increasing functions of », and cannot, for 
ppX<2a+f8, exceed a fixed value Ny obtained by evaluating them for 
p=2a+f, ie. «-=0. Thus, when N>N, we must have l/e ~ N—N), 
which is of the order of magnitude of NV, a large number. But, when 
N<N,, (5) can be satisfied by taking «=c(1—N/N,)?, i.e. of the order 
unity. This is because the leading term in the first sum then behaves 
like a ‘ generalized Zeta Function’ as in the usual theory of the B.E. 
condensation, and can be expanded in a power series in ¢«!/?, of the 
general type Nj(l—ae!/?-+...) (Robinson 1951). In other words, as 
the density increases (at constant temperature) a point is reached at 
which the chemical potential ceases to rise visibly, but becomes almost 
constant. This represents a condensation process; physically it means 
that all subsequent particles beyond Ny, go into the lowest energy level. 
Since we are always keeping N finite, there is no question of a mathematical 
discontinuity in the gradient of 1 say, but we have an abrupt change, 
over a small interval, which would be physically indistinguishable from 
a true singularity. 

To find the conditions for the condensation to occur, we must solve 
the equation for Np, 1.e. 


[2B T 2 (/3 ) | 7/20 «( /e 2 
= — j(2a+ B) se | tee Ih ee eZ) > ay 
of ved we JB 4 jn/op+1 J& 
War (i) 


where the right-hand side is a function of temperature through « and p. 
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The exact evaluation of the various sums in (7) is exceedingly complicated. 
However, if we write it in the form 


SNoxBe™ 7 a jt2a+A) Oe Ibs (/5-/5) 
a j us 4 iP 
{2x . i| 2 
iff a om( [eff ran oe by) 
TT j=n/2B+1 y 


the left-hand side is a function only of the number density of particles, 
N,/L2D, and of the temperature, whilst the right-hand side may be 
expanded as a function of « and £, which we take to be small. An 
outline of the method of expansion is given in the appendix. The 
result may be written 


(7',/T)32=1+B1/2(—0-19+ 0-648 In B—0-432 In «) 40-4320"? In B, (9) 
where (NLA D\hF/2aM kT oA = 26125 os ee ee 


ie. T', is the temperature of the B.E. condensation in an infinite cube, 
at density N,/L?D, according to the usual theory. 

This formula may be expressed in a more self-evident manner if we 
use the spacing of atomic layers in the film as our unit of length. This 
unit is 31/42-1/3(L2D/N,)/% (Frederikse and Gorter 1950), which is 3-72 A 
in ordinary liquid He. Then we have «=1-17(7,/T)L-?, B=1-17(T,,/T)D ~~, 
i.e. «1? and B-!/? are approximately the length and thickness of the 
box in numbers of atomic layers, when the temperature is near the 
lambda point of the bulk fluid. In these units we have 


T\ 3/2 T/T ie T 
(7) -1-3(7) {0-15-o-504 In #) 1-912 In D—1-008 In Z} 


1-008 /7',\ 1/2 
-" (Fz) In D. cn. cite. cel Oh ee 


§3. SHAPE-DEPENDENCE OF CONDENSATION TEMPERATURE 


Equation (11) gives the temperature, 7', of the B.E. condensation in 
a block of dimensions L x L x D in ‘ layer units ’, in terms of the transition 
temperature 7’, in an infinite cube of bulk liquid of the same density. 
Now let us consider some special cases. 

(a) If L > «, whilst D remains constant, it is obvious that 7’ tends 
to zero. This result is not a spurious consequence of the approximations ; 
it can be seen to be due to the second sum in (8) which has a singularity 
at a0. In the limit of large NV, it is easy to show that 


(T/T)? x 1/B1/2 In N. 


This corresponds to Osborne’s result for a two-dimensional gas, but now 
we see that every layer of finite thickness shows the same behaviour, 
the difference being that the thicker the system the larger must be the 
value of V to achieve a given lowering of the condensation temperature. 
Osborne's critical thickness occurs at B=1, when the system does indeed 
have exactly the properties of a two-dimensional gas, but in the limit, 
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as N + oo, this is irrelevant, and all layers of finite thickness behave 
alike. : 

(6) If both L and PD tend to infinity together, in any way, 7'> 7, 
and we have the usual theory of the bulk fluid. 


(c) If L and D are equal and finite, two terms in (11) cancel. The 
resulting formula is not quite correct, however, because we have assumed 
throughout that «<B, ie. D<L. The necessary adjustments are éasily 
made, and we obtain 

(TL)? = 1—L 47/1) {— 1-15—0-504 In (1/7) 4+-1-512 nL}. (12) 

This formula is, of course, an approximation valid only for large L; 
we need only consider the case In L>1-15/1-51, where it is evident that 
T>T,. The maximum possible value of 7/7, is 1-25, and the ratio 
tends rapidly to unity as L increases. Thus, the lambda point in a 
finite cube is slightly higher than in the infinite assembly, except when 
the sides of the cube are only a few angstréms in length, i.e. when it 
contains only a few score atoms. 

(d) If D<L, and L is constant, the dominant factor on the right-hand 
side of (11) is that containing In L, and it is clear that 7/7’, decreases as 
PD decreases. The trend is monotonic, and 7 +0 as L—0. (One might 
object that we have then violated the condition that B be small, used 
in evaluating (9). Study of (8) for the case B>7/2 shows that the first 
sum may be neglected, and the second sum behaves like fh! Ina, the 
same dominant term.) 


$4. THE ‘ SCATTERING LENGTH ’ 


It is evident from the above analysis that we can make no progress 
in studying the lambda-point phenomena in real films of helium if we 
insist that they shall form parts of infinite assemblies. If we were to 
consider every observed helium film as a complete assembly, we should 
find the condensation temperature varying with the size of, say, the glass 
plate on which it is formed, which is not in agreement with observation. 
In any case, the calculated values of 7 from formula (11) turn out to be 
much smaller than the observed values, if we put in the macroscopic 
lengths of the apparatus, of the order of | cm, for L. 

The only way out of this difficulty is, apparently, to make the following 
new assumption: a macroscopic assembly of liquid helium must be 
treated as if it were made up of independent sub-assemblies, none of 
whose linear dimensions may be greater than a fixed length, L). Thus, 
for the bulk liquid we consider a cube, of side Ly, with the particles 
confined strictly within its boundaries, and construct the macroscopic 
volume by piling up a number of identical, physically independent 
cubes. When a film is thinner than Lo, we take rectangular blocks, of 
side Ly, to cover the area. e 

The effect of this assumption on the theoretical properties of the bulk 
liquid is negligible, so long as L, is moderately large. We then have 
case (c) above, and it is clear that the effect is merely to raise the predicted 
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lambda point by a small amount. Since the theoretical condensation 
temperature of a B.E. gas of the density of liquid helium is, in any case, 
about 50% larger than the observed value, this further small change is 
quite unimportant. 

On the other hand, by choosing a suitable value of Lo, we may be able 
to fit the observed variation of the condensation temperature with film 
thickness. That is, we solve (11) as in case (d) above, and interpret D 
as the thickness of the film. 7 cannot be taken at its face value as the 
transition temperature in the film, because of the discrepancy between 
the observed and calculated values of 7’, for the bulk fluid. It is reasonable 


0-8 


0-7 


06 
4 6 8 10 


D (Atomic Layers). 
Transition temperature as function of film thickness. 


observed —-— — calculated 


to assume however, that the cause of the discrepancy acts proportionately 
in bulk fluid and in film, and we can ‘ renormalize ’, by comparing observed - 
and theoretical values of 7/7',. If we insert the experimental values 
of D and 7/7, (Brewer and Mendelssohn 1953) into (11) we obtain a 
series of values of L, ranging between 150 and 250 atomic layers. This 
represents a variation of only 10° in In ZL, which is the significant 
parameter in the theory, whilst D goes from 4 to 10 layers. The value 
Ly= 200 layers=7-4 x 10-& em has been chosen as the best estimate for 
the scattering length, and in the figure the theoretical variation of the 
transition temperature with the thickness of the film is compared with 
the experimental curve. 
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It will be seen that the agreement is remarkably close, particularly 
in the central regions of the curves. Ly is the only adjustable parameter 
in the theory, yet the curves can be made nearly to coincide over quite 
a large range of values. We should not be surprised to find discrepancies 
at the lower end, because there the theory is not properly valid, terms in 
B ~ D™ having been neglected. 


§ 5. ConcLUSION 

It is clear that by introducing a standard ‘ scattering length’ of the 
order of 7 10-8 cm into the theory of liquid helium we can improve the 
usual theory in two entirely distinct ways. We avoid the paradox that 
a physically thick layer, of infinite extent, would behave differently from 
a sample of an infinite cube and we obtain a very good theoretical estimate 
of the variation of the transition temperature with film thickness. 

There are other considerations which lend some support to this 
hypothesis. The scattering length is of the same order of magnitude as 
the thickness of the saturated helium film below the lambda point, i.e. 
about 3x 10°%cem. Now we have shown (eqn. (6)) that the chemical 
potential of the helium atoms below the lambda point is given by 
kT'(2a+8). «is fixed by our choice of Ly, but 8 depends on the thickness 
D of the film, and has a minimum value, «, when D=JL,. Thus, there is 
a distinct gain in negative free energy if the film thickens to this extent— 
but no further gain from any greater thickness, since we should then have 
D>L,, which contradicts our assumptions (or, rather, we should have to 
construct new fundamental cubes on the surface of the layer). This 
minimum chemical potential is the same as that in the bulk liquid, as 
it must be in a saturated film. Above the lambda point, the small gains 
obtained in the chemical potential in this way are outweighed by large 
entropy terms, and there is no reason to depart from the normal theory 
of adsorbed gases, which gives a thin film. 

Another difficulty of the standard theory which is avoided concerns 
the interpretation of the ‘ condensed’ atoms as superfluid atoms, for if 
these are assumed to lie in the lowest energy state of an infinite container, 
they would have zero momentum, and would be incapable of flow 
(Mendelssohn 1945). But, in our interpretation, the ‘ condensed ’ 
atoms lie in the ground state of a box of side LZ, and it can easily be 
shown that this corresponds to their having a ‘ zero-point velocity ’ of 
70 cm/sec, if Ly is about 7x10-®cm. This is of the same order of 
magnitude as the observed maximum velocity (30-50 cm/sec, Bowers 
and Mendelssohn 1952), for flow through superleaks, and provides a 
simple explanation of the fact that this velocity does not seem to depend 
on the dimensions of the apparatus. 

The agreement may well, of course, be fortuitous. Leaving aside the 
question, to which we shall return below, whether this equating of 
zero-point velocity to superflow velocity is tenable, we are still no nearer 
to explaining why the flow of the condensed atoms should be frictionless. 
To discuss this, we should need to allow exchange of particles between 
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neighbouring boxes, and it is difficult to see how this can be formulated 
in a consistent manner. But in this respect we are no worse off than in 
the usual B.E. gas theory; all that we have done is to assign a new 
distribution of energy levels to the whole assembly by insisting that the 
wave functions shall vanish at each of a set of plane boundaries. 

It has recently been suggested (Brewer and Mendelssohn 1953) that 
the anomalous behaviour of unsaturated helium films below the lambda 
point of the bulk fluid may be explained by the formation of ‘ clusters ; 
at least 10-&cm in diameter, of the bulk fluid. This is supported, at 
least qualitatively, by the present theory, since we provide a reason for 
the clusters to grow to a certain thickness Ly, but not larger. 

Throughout this paper, the object has been to show that the usual 
theory of B.E. condensation is paradoxical, and cannot be applied consist- 
ently to the case of thin layers. We have shown what modification is 
necessary to resolve the paradox and to give some agreement with the 
observed properties of helium films. It is worth discussing briefly 
whether there is any a priori justification for the assumption we have 
made—the introduction of a standard length. 

In treating liquid helium as an ideal gas, we have neglected two 
important factors: the van der Waals attractive forces between the 
atoms and their finite impenetrable volumes. The attractive forces, 
although responsible for the condensation to the liquid state, are very 
weak, and it seems reasonable to take account of them by assuming that 
inside the liquid they give rise to a smeared-out uniform potential, which 
rises sharply at the boundaries and reflects the atoms back into the 
occupied volume. This model gives a good theoretical estimate of the 
vapour-pressure curves of ?He and *He (Daunt 1952). But it cannot be 
correct to treat the atoms as if they were free to bounce unchecked from 
side to side of the whole container. It is true that the density of liquid 
helium is low, and the spacing between the atoms much larger than the 
gas-kinetic diameter, but there must be collisions and scattering of one 
atom by another. After having gone a certain distance, an atom will 
have bs forgotten’ its original momentum, owing to this scattering. 
That is to say, it cannot be correct to discuss the energy levels as if they 
were those of the whole container, the only conditions on the wave 
functions being that they should vanish on the walls. The picture of 
plane waves persisting unchanged from one side of the box to the other 
must be modified to allow for this incoherent scattering. Our assumption 
is tantamount to simulating the scattering by inserting walls, distant 
Dc; apart, and insisting that the wave functions shall vanish on these walls 
also. The detailed shape of these elementary volumes is not important ; 
the asymptotic distribution of the eigenvalues of the wave equation 
depends only on the general dimensions of the assembly (Courant and 
Hilbert 1931, p. 373) 

The distance Ly) may thus be interpreted as a ‘scattering length ’. 
But if we attempt to calculate it directly, from kinetic theory, it turns 
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out to be far smaller than the value we have adopted. It will, of course, 
be rather larger than the mean free path for collisions, as calculated from 
the geometric cross-section of the atoms, owing to the persistence of 
velocity ; the effect of an impenetrable wall on an impinging particle 
is much more drastic than that of a glancing collision with a like particle, 
so that the walls must be rather further apart to give the same scattering 
effect. But even then, about 10-7 cm would seem to be an upper limit, 
at the density of liquid helium. 

That such a calculation must be fundamentally correct is shown by 
the fact that the ‘ zero-point energy’ in liquid helium is large, corres- 
ponding to a velocity of about 104 cm/sec, whereas, on our model, as 
pointed out above, this would only be about 70 cm/sec. This is a 
fundamental difficulty in all theories based on the B.E. condensation 
(Mendelssohn 1945). Since the zero-point motion is just ‘ uncertainty 
momentum ’, the velocity is inversely proportional to the scattering 
length, for which a value of about 5x 10-* em is thus suggested. 

Thus, although we have taken a step in the right direction to bring the 
B.E. gas model nearer to the physical reality of liquid helium, there is 
still a large gap which cannot be bridged ; if L, is of the order of a few 
angstréms, the condensation phenomena are inexplicable—if it is as large 
as 7x10-8em, the model has no firm basis in kinetic theory. We 
apparently need a theory in which lengths of both these magnitudes 
occur, in order to explain, on an equal footing, all the phenomena. 


IT am indebted to Dr. K. Mendelssohn and Mr. D. F. Brewer for several 
valuable discussions, and for the experimental data shown in the figure. 


ARE BeliN el) TeX. 
The evaluation of the various sums occurring in (8) is carried out 
along the following lines. A typical term would contain, say, 
n/2B ‘ 
U2) 9-*e”, 
‘ j=l 
where b<f. Thus, jb<7/2 for the values of j included in the sum, so 
that the exponential term may be expanded in a power series. We have 
n/2B Re . ° fi 
UAL Lyr(joys/s! 
j=1 s=0 


(ve) n/2B : 
BOSSA Zigee 
s=0 j=l 


But, by the Euler-McLaurin theorem (Whittaker and Watson, p. 128), 


30/28 7/28 
os fe SS: | 


j=1 


ee) (—1)=3 (2t-1 71/28 
aot da +4 (1+ (n)26)° V+ 2 SB, Ee ‘| . 
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We need only retain the lowest power of 8, which comes from the 
upper limit of the integral, i.e. 


E {b§/(s—r+ 1)s !}(a/2B)s-* +2 


U,~ 
8 


= (2B/m)"-1 3 (mb/28)%/(s—r+ ls}, 
s=0 


and the series converges rapidly. 
For the leading term of (8), for which r= 3/2, it is better to use 


7/28 oa) 

Ej a8) Bi, 

j=1 j=n/2B+1 
giving a correction term, of order 61 to the known Zeta Function. 
The logarithmic terms arise when s—r=—1 in the above expansions and 


the series can then be evaluated by the theorem of Euler’s limit. 

In principle, the calculation can be carried to higher orders of approxi- 
mation, using the asymptotic Kuler-McLaurin expansion, but it becomes 
exceedingly complicated. The most important terms in (9) are those 
involving Ing and Inf, and these have coefficients which are exactly 
2/4/76(3/2) and 3/./7¢(3/2) respectively. The number 0-19 is the resultant 
of many sums and differences of terms like U,, but it evidently would not 
affect the theory very greatly if it were a few per cent in error. 
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LIX. CORRESPONDENCE 


Anomalous Surface Tension in Heliwm II 


By D. F. Brewer and K. Menpexssonn, F.R.S. 
The Clarendon Laboratory, Oxford 


[Received March 17, 1953] 


In trying to interpret recent observations on the formation of the 
helium IT film (Brewer and Mendelssohn 1953), we have suggested the 
presence of bulk liquid at pressures well below the saturation pressure. 
In order to account for such a phenomenon, a quantity of the dimension 
of a length, d, characteristic for liquid helium IT, was introduced. The 
significance of d is that of the de Broglie wavelength of the lowest energy 
state, i.e. the distance within which the position of a helium ‘atom in 
this state cannot be determined. It may then be energetically 
advantageous for clusters of bulk liquid to be formed of the order of 
magnitude of d. It was also concluded from the experiments that d 
would be at least ~10-®cm. The existence of such a characteristic 
mean free path in liquid helium II has since been made plausible on 
theoretical grounds by Ziman (1953). We have further concluded that 
such a model should lead to an anomalous surface tension effect. The 
_ macroscopic surface tension does not show a discontinuity at the lambda 
point (Allen and Misener 1938), as is indeed to be expected owing to the 
absence of a latent heat. However, if our model is correct, one can expect 
the surface tension to become size-dependent when dimensions of the 
order of d are approached, because of the extra energy available for 
liquid clusters of this size. An experiment was therefore designed to 
investigate this question. 

A cylindrical container C (fig. 1) 12 cm long and of internal diameter 
2-3em was closely packed with about 1000 plane circular glass plates 
with an estimated gap of ~10~4 cm or less between neighbouring plates. | 
At the bottom of C a glass tube 7’ was attached. The whole container, 
which could be filled with gas through the line L, was immersed in liquid 
helium. 

Helium gas was introduced into C just below the lambda point until 
some bulk liquid appeared in 7’. The temperature of the helium bath 
was then raised to just above the lambda point. In doing this, some 
liquid should evaporate from 7’, but as the gas space as well as the 
relative increase in temperature were small, this effect was insignificant. 
In addition, one. must expect the helium film of ~10~* cm thickness 
to disappear from the surface of the glass plates, and this should lead to 
a small rise of the liquid level in 7’. 
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An increase of bulk liquid in the tube 7’ was indeed observed, but 
the amount exceeded that calculated for the film by between 10 and 
20 times. 

Since ordinary capillary condensation due to the normal surface 
tension will change very little in the small temperature interval (<0-1°) 
through which the temperature of the container C was raised, the 
appearance of this large amount of bulk liquid must be connected with 
the lambda phenomenon. It seems indeed that below the lambda point, 
a large amount of liquid helium IT was held in the gaps between the 
plates. As the lambda temperature is exceeded, the significance of d is 
lost, and with it the extra energy holding the liquid in clusters of small 
dimensions. 


While thus our experiments indicate that the postulated effect exists 
certain features of it require further elucidation. For instance another 
appreciable quantity of bulk liquid appears as the temperature is raised 
from 2-2° to 2-5°x. It is not clear whether this is due to the fact that 
at the lambda point, further expulsion of liquid from the plates becomes 
very slow, or whether the effect actually extends into the region of the 
slow anomalous increase of entropy between the lambda point and 2-6°K. 
Moreover, it may be significant that in the present as well as in the 
earlier experiments, the plane of the solid surface on which the helium 
was deposited was horizontal. Finally, the relation between the 
clustering effect and the well known film phenomenon is not quite clear. 
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The most likely explanation seems to be that the tendency to form 
liquid clusters of the order of the de Broglie wavelength is the primary 
effect. of which the film presents only one particular aspect 
‘The existence of an anomalous surface tension effect in liquid helium IT 
will go far in explaining a number of curious features in earlier experiments, 
It has been shown (Bowers and Mendelssohn 1950) that the anomalously 
high transfer rates sometimes observed (de Haas and van den Berg 1949) 
were due to surface contamination by solidified gases, and that surface 
roughness in general leads to higher transfer rates (Chandrasekhar and 
Mendelssohn 1952). It is now clear that cracks in surfaces in contact 
with liquid helium II are filled with bulk liquid. The new effect. is 
probably also responsible for the unexplained high adsorption of helium 
below the lambda point observed by Long and Meyer (1949) in packed 
rouge powder. The peculiar transport phenomena observed by the 
same authors (1952, Method II) in helium films of only 14 statistical 
layers may also possibly be due to anomalous capillary condensation in 
the superleak over the surface of which the flow was thought to occur. 
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Flow of an Inviscid Supersonic Jet past an Infinite Wedge 


‘Bye DG Pack and L. Roperts 


Department of Mathematics, University of Manchester 


{Received 18 March, 1953] 


THe general character of the two-dimensional inviscid flow, past an 
infinite wedge, of an infinite stream with supersonic velocity at infinity 
in the upstream direction, is well-known. When the wedge is too blunt 
or the stream not sufficiently fast, there is a shock-wave astride the 
~ flow at some distance upstream from the wedge. On the line of symmetry, 
between the central point of the shock-wave and the tip of the wedge, 
the speed is subsonic, and at the tip itself there is a stagnation point, 
When the wedge is narrow enough or the stream sufficiently fast, the 
shock-wave system is attached to the tip of the wedge. There is a small 
range of wedge angles and Mach numbers for which the flow between the 
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shock-wave and the wedge is subsonic, but almost always when the 
shock-wave is attached the flow up the sides of the wedge is uniform 
and supersonic and the shock-wave straight. 

Although there have been many papers on the subject of flow past 
a wedge with a finite length of sloping edge in an infinite stream, the 
authors are not aware of any discussion of the flow, past an infinite wedge, 
of a stream of finite thickness such as a jet. Examination shows that 
this problem has some interesting features. 


UNIFORM SUPERSONIC 


FLOW FROM 
LOW FROM NOZZLE PLANE OF SYMMETRY 


Consider (see the figure) an infinite wedge impinged upon by a uniform 
jet of gas flowing with supersonic speed (Mach number /,) at atmospheric 
pressure from infinity. The speed is supposed to be great enough to 
produce an attached straight shock-wave with supersonic flow up the 
sides of the wedge if placed in an infinite stream with the same speed. 
The shock-wave will be straight for the jet also; for, upstream of the 
characteristic AB from the point at which the shock-wave meets the 
boundary, the fact that the jet is not of infinite width can have no 
influence. The shock-wave deflects the bounding stream-lines and 
increases the pressure. This increase must instantaneously be removed 
by means of a Prandtl-Meyer expansion down to atmospheric pressure, 
This is indicated in the upper half of the figure, by the pencil of Mach 
lines from A (the broken Mach lines denote expansions, solid Mach lines 
compressions). 

It will be observed that the situation behind the shock-wave is exactly 
that which occurs when a jet emerges with uniform supersonic velocity 
from an orifice, the velocity in this case being the velocity (corresponding 
to Mach number M,) behind the attached shock-wave. The flow is 
uniform up to AB. The subsequent motions on the upper and lower 
faces of the wedge are respectively the upper and lower halves of the 
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supersonic jet described above, and are therefore quasiperiodic with 
shock-waves occurring in a well-known manner (Pack 1948). (The 
figure illustrates one complete period of the jet structure, assuming 
absolute periodicity. ) 

Professor Lighthill has pointed out to the writers that an experiment 
designed to reproduce the state of affairs outlined above would be of 
great interest. The orthodox supersonic jet is quite well understood, 
although not many experiments have been carried out on the two- 
dimensional form. The flow up the wedge side considered in this note 
contains a half-jet bounded, not by another half-jet, but by a solid wall. 
The usual plane of symmetry in the gas has been replaced by the solid 
boundary of the wedge. This solid boundary introduces a boundary- 
layer effect, and it would be profitable to see what differences are induced 
in the flow by its presence, by comparing studies of the flow in this 
situation with jet flows obtained by the usual means. 

Theoretical work is being undertaken on other and more difficult 
problems relating to the flow of supersonic jets past bodies when the 
bow shock-wave is detached. 
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ERRATUM 


The Role of Preferred Orientation in Elasticity Investigations, by G. 
BRADFIELD and H. Pursey, 1953, Phil. Mag., 44, 439. 
Equation (2) should read : 
— 125¢(¢,,+2¢,2)(K —1/5) 
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LX. Notices of New Books and Periodicals received 


Thermal Diffusion in Gases. By K. E. GRew and T. L. Less. [Pp. xi + 143.] 
(Cambridge Monographs on Physics ; C.U.P. 1952.) Price 22s. 6d. net. 
Tue need for a book of this character has been evident for some time. Most 
of the large and scattered body of experimental and theoretical work on the 
thermal diffusion effect in gases and liquids, up to the end of 1951, is here 
mentioned and correlated. The study of thermal diffusion as a means of 
investigating intermolecular forces is most fully dealt with and the book 
will become the standard monograph on this study, alongside de Groot’s 
monograph on the Soret Effect and the Furry and Jones (1946) paper on the 

Clusius—Dickel column. , 

The first two chapters are devoted to an outline of the discovery, nature 
and theoretical treatments of the thermal diffusion effect in gases. Relevant 
parts of the Chapman—Cowling theory are presented in detail and developed 
in terms of various molecular models. The third and fourth chapters cover 
the experimental methods and results of ‘two-bulb’ experiments, and in the 
fifth chapter an excellent and detailed discussion and comparison is made 
of these results with the theory. An account of the diffusion-thermoeffect is 
given in the next chapter. This is most welcome, as it is probably the first 
detailed description of Waldmann’s work to appear in English, and as it also 
includes a phenomenological theory of the effect. The Clusius—Dickel column 
is dealt with next, necessarily rather briefly, in 16 pages. The simple theory 
of the operation of the column is given, together with some examples of its 
use for concentrating isotopes. It is a pity that no reference is made to Clusius’ 
production of 99-8°9,N, nor to Jensen’s theory of the column which is that 
usually quoted by Clusius. The last chapter covers the Soret Effect, again 
briefly, and includes work published since the appearance of de Groot’s mono- 
graph in 1945. The appendices include a useful list of values of « for nearly all 
mixtures examined experimentally. R.H.D. 


The Method of Descartes. (A study of the Regulae.) By L.J. Brox. (Oxford : 
Clarendon Press, 1952.) [Pp. 316.] Price 30s. 

THis work is a thorough commentary on the Regulae and to a great extent on. 
the Discours. Mr. Beck argues, among other things, that the latter cannot 
properly be understood without the former. His exposition of Descartes’ 
ideas of method and of their importance for Descartes’ works in philosophy, 
mathematics and physics is lucid and illuminating. Mr. Beck’s book will, 
the reviewer believes, prove to be of real use to those who wish to understand 
the guiding conceptions of Descartes. 8. K. 


British Association Mathematical Tables, Vol. X, Part 2.—Bessel Functions of 
Positive Integer Order. By W. G. Bicktny, L. J. Comriz, J. C. P. Minter, 
D.H.Sapier and A. J.THompson. [Pp. xl+255.] (Cambridge University 
Press, 1952.) Price 60s. ; 

THE publication of this volume is an event of major importance to table-users. 

Vol. X, Part 1, containing Bessel functions of orders zero and unity, appeared 

in 1937; Part 2, now issued, contains the functions J,, Y,, J, and K, for 

further integral n up to 20. The delay in fulfilling a scheme conceived so long 
ago is easily understood, not only by reason of the size of the task but because 

of many unavoidable interruptions. However, the authors have produced a 

set of tables of lasting value, and which in design and presentation reach the 

highest standard. R. AaB 


[The Editors do not hold themselves responsible Jor the views 
expressed by their correspondents. | 
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